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1.      INTRODUCTION 

The  purpose  of  this  paper  is  to  investigate  and  recommend  suitable 
types  of  measurements  on  the  lunar  surface  which  will  result  in  the 
establishment  of  a  selenodetic  control  network.      In  this  paper,   we  will 
define  a  selenocentric  celestial  coordinate  system,    in  which  we  may 
make  measurements,    and  in  turn  determine  an  astronomic  position  on 
the  lunar  surface.      Combining   this  astronomic  position  with  a  knowledge 
of  the  deflection  of  the  vertical  at  this  position,   we  could  then  determine 
the  coordinate  of  this  position  in  a  selenographic  system. 

Therefore  this  paper  will  be  divided  into  two  primary  areas  which 
are: 

(1)  Establishment  of  a  selenocentric  celestial  coordinate    system 
and  measurements  in  that  system  from  the  lunar  surface. 

(2)  Determination  of  the  deflections  of  the  vertical  on  the  lunar 
surface. 


2.      ESTABLISHMENT   OF  AND  MEASUREMENTS  FN  THE    LUNAR 
ASTRONOMICAL  COORDINATE   SYSTEMS 

On  the  earth  most  of  the  geodetic  astronomic  calculations  are  made 
using  a  geoequatorial  system  where  the  star  coordinates  are  right  ascen- 
sion and  declination.      There  are  a  great  number  of  catalogs,   maps  and 
almanacs  giving  the  coordinates  of  the  stars  and  other  celestial  bodies 
in  this  system.      If  the  coordinates  in  this  system  could  be  easily  trans- 
formed to  a  lunar  equatorial  system,   all  the  methods  of  geodetic  astronomy 
could  be  emploj^ed  on  the  moon  for  selenodetic  purposes.      In  this  section 
we  will  show  such  a  transformation  complete  with  error  analysis,   while 
also  discussing  such  related  factors  as  time,   location  of  lunar  pole  stars 
and  observations  in  the  horizon  system. 

We  will  also  show  some  methods  which  may  be  used  on  the  lunar 
surface  to  determine  astronomic  position.      These  methods  will  be  dis- 
cussed in  relation  to  the  instrument  requirements,    and  accuracy  of  the 
system  as  a  whole.      Included  among  these  methods  is  an  inertial  system, 
which  is  presently  installed  on  the  lunar  module,   and  which  in  itself  can 
establish  a  position  on  the  lunar  surface. 

The  methods  and  instruments  presented  are  not  assumed  to  be  a 
complete  listing,  but  are  based  on  practical  factors  which  may  obviously 
limit  other  methods  because  of  physical  size  or  complexity  of  equipment, 
time  required  to  achieve  an  astronomic  position,   or  excessive  burden 
being  placed  on  the  astronauts,    either  during  training  or  during  the 
operation  itself.      In  this  light  we  would  most  like  to  find  a  method  which 
is  as  conservative  of  the  astronauts  time  as  possible  (in  both  training  and 
employment),   achieves  satisfactory  results,   and  does  not  involve  undue 
calculations  or  recording  procedures.      In  conjunction,   we  wish  the 
instruments  involved  in  this  method  to  be  simple,   portable,   and  of 
course  accurate. 


2.1     Lunar  Astronomy 
Kolaczek  [14]  defines  a  selenoequatorial  coordinate  system   (lunar 
right  ascension  system)  where  the  axis  of  rotation  of  the  moon  defines 
the  celestial  pole,   which  is  the  basic  direction,   and  the  celestial 
equator,   which  is  a  plane  perpendicular  to  the  rotation  axis,   containing 
the  center.      The  latter  is  defined  as  the  basic  plane.      This  system  is 
affected  by  precession  and  nutation  (physical  libration)  of  the  lunar  axis 
of  rotation.      The  selenoequatorial  system  is  analogous  to  the  geoequatorial 
system  on  earth,  but  the  selenoequatorial  coordinates  of  the  stars  change 
their  values  faster  than  the  geoequatorial  coordinates  due  to  the  lunar 
precession  being  1,360  times  faster  than  the  earth's.      The  selenoequatorial 
coordinates  are: 

oc       -     Lunar  right  ascension  measured  on  the  lunar  celestial 
equator  from  its  ascending  node  on  the  ecliptic.      It  is 

o  o 

measured  to  the  east  from  0    to  360  . 
d       -     Lunar  declination  which  is  the  angular  distance  from  the 

lunar  celestial  equator.      It  is  measured  from  0    to  90  , 

positive  in  the  northern  celestial  hemisphere,   and  negative 

in  the  southern. 
Now  we  will  take  a  coordinate  system  with  known  coordinates  and 
transform  these  coordinates  into  the  selenoequatorial  coordinate  system. 
The  translation  of  any  coordinate  system  from  one  point  in  space  to 
another  changes  the  values  of  the  spherical  coordinates  of  fixed  points 
on  the  celestial  sphere.      This  change  is  caused  by  the  parallax  or  the 
translation  of  the  origin  of  a  coordinate  system  and  the  aberration  caused 
by  the  different  motion  of  this  newly  translated  system.      The  translation 
of  mean  earth  oriented  coordinates  into  mean  lunar  oriented  coordinates 
does  not  require  consideration  of  parallax  if  we  consider  only  the  coordi- 
nates of  stars,   as  we  may  omit  the  parallax  of  the  earth-moon  distance, 
due  to  the  great  distance  of  the  stars,   as  compared  to  the  earth-moon 
distance.      The  transformation  of  the  apparent  coordinates  however  requires 
consideration  of  the  influence  of  the  moon's  aberrations. 
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2.1.1  Transformation  of  Mean  Geoequatorial  into  Mean  Selenoequatorial 
Coordinates 

The  transformation  of  mean  geoequatorial  into  mean  selenoequatorial 
coordinates  can  be  performed  in  two  steps  [14].      First  the  mean  geoequa- 
torial coordinates  ex. ,    5  are  transformed  into  the  mean  ecliptic  coordinates 
X,   j3  by  the  well  known  formulas: 

sin#  =  cos  6  sin6  -  sine  cos 6  sintt 

cos/3   cosX  =  cos  6  cos  Q! 

cos/3  sinX  =  sine  sin 6  +  cos  e  cos 6  sin<2 
where  €    the  obliquity  of  the  ecliptic  is 

€  =  23°.  452294  -   0°.  0130125  T  -   0°.  00000164  T2   +  0°.  000000503  T3 
where  T  denotes  the  time  measured  in  Julian  centuries  of  36,525  ephemeris 
days  from  the  epoch  1900  January  0.5ET  =  J.  D.    2415020.0. 

Next  the  transformation  formulas  of  mean  ecliptic  coordinates  into  mean 
selenoequatorial  coordinates. 

sin  dm     =  sin/3  cos  I  +  cos/3  sin  I  sin  (  X  -  £1) 

coscL  sin  ex  n  =  sini3  sin  I  -  cos/3  cos  I  sin  (\  -  Q)    f  2.1 

cos  dm   cos  a  a  =  -  cos/3  cos  (  X  -  &) 
where  I  is  inclination  of  the  lunar  equator  to  the  ecliptic 

I  =  1°  32'  01" 

d     The  mean  longitude  of  the  ascending  node  of  the  lunar  orbit  on 
the  ecliptic  is 

Q=  259°.  18305556   -  1934°.  14200833  T   +  0°.  002078  T3   +  0°.  000002  T3 
where  T  is  defined  as  above. 

2.1.2  Transformation  of  Mean  Selenoequatorial  Into  Apparent  Selenoequatorial 
Coordinates 

The  transformation  of  mean  into  apparent  coordinates  is  given  in 
a  form  [14]  similar  to  that  used  for  the  earth's  equatorial  coordinates 

a  sapp     =    Of  m     +  Ai  ai  +  Bi   bi  +  Ci  cx  +  Di  dx  +  ra  Ua 

d  sapp     =     d  m  +  Ai  al  +  Bi  bl  +  Ci  ci  +  Di  dl  +  Ta  Ud 
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Similarly  to  Bessel's  equations  on  the  earth,   Ai  and  Bi  are  the  effects 
of  precession  and  nutation  to  date  from  the  determination  of  the  mean 
coordinates.  The  Ci  and  Di  terms  are  the  effect  of  lunar  annual  abberation. 
The ra  u  /  d  /a  \  terms  are  the  corrections  for  proper  motion.        The 
ai  ,bi  ,ci,di,ai  ,bl  ,cl  ,dl  coefficients  are  the  lunar  Besselian  star  constants. 
Specifically: 

Ai  =  -  (Pot  +    cr  )sinl 

Bi  =     p 

Ci  =  -  K  sin  (Lap  -  d+  180°)  cos  I 

Di  =  -   K  cos  (Lap  -  Q  +  180°) 


where: 


Po" 


p-.ff 


K,  k 


Ls  ,  ll    ,     La  ; 


V 


Daily  motion  of  the  ascending  node 
of  the  moon's  orbit  on  the  ecliptic 
(  =  -  0°.  0529539222) 

-  Number  of  ephermeris  days  from  the 
epoch  to  of  the  mean  coordinates  &   n  ,  dB 

-  Inclination  of  the  moon's  equator  to  the 
ecliptic   (  =  1°  32'  01"  ) 

-  Physical  librations  in  inclination  and  in 
ascending  node  respectively  (see  Appendix  #3) 

-  Constant  of  the  earth's  and  of  the  lunar 
aberration  respectively 

K  =  k(l  -  Vi  cos(Ls  -  li)/Ve) 

(k  =  20". 4958) 

Longitude  of  the  sun,   the  moon  and  the 

moon's  apex  in  the  motion  around  the 

sun  respectively 

LaP  =  Ls  -  90°  +  A  A 


A  A  .      "I    fVl  siD    (  L.   -    ll  )  \ 

AA  =  ta°      U.-VicobCL.-!!); 


Ve,Vi  -     Velocity  of  the  earth's  and  the  moon's 

orbital  motion,    respectively 

(Ve  =  29.8  km/sec 
Vi   =1.02  km/sec) 

The  lunar  Besselian  star  constants  are: 
ai  =  (  cot  I  -  sin  ol  B  tan  d0  ) 

,  a 

ai  =  -  cos  ol  B 

bi  =  cos  Q!  B  tan  dm 

bi  =  -  sin  ol  n 

ci  =  -  cos  Q!  B  sec  dB 

t  0 

ci  =  sin  ol  n   sin  dB  +  tan  I  cos  dB 

di  =  sin  a  B  sec  dB 
di  =  cos  a  B  sin  dB 

The  proper  motion  will  take  on  similar  values  to  those  computed 
on  the  earth: 

ra  =    t 

/4,  Md    =     The  components  of  proper  motion  of  the  star 

2.1.3     Lunar  Horizon  System 

Measurements  on  the  lunar  surface  will  most  likely  be  made  in  the 
horizon  system.    In  the  moon's  horizon  system,    altitude  aih  and  azimuth  Aih 
are  defined  in  the  same  way  as  that  of  the  earth's  horizon  system.      The 
transformation  between  the  lunar  horizon  system  and  the  lunar  selenoequatorial 
system  similar  to  that  given  by  Mueller  [16]  is  as  follows: 

cos  d  sin  (  LHA~-.A-  ol      )  =  -  sin  Zi  sin  Aih  ^ 

sin   d  =  cos  Zi  sin<£  +  sin  Zi  cos  A  cos^  >2.2 

cos  d  cos  (LHA^  -A-  a      )  -■=  cos  Zi  cos£-  sin  Zi  cos  Aih  sin<£J 


where  LHA0  -  The  lunar  hour  angle  of  the  node  measured  from  the 

La 

reference  median 
<3?  -  The  lunar  astronomical  latitude 
X    -  The  lunar  astronomical  longitude 

Zi  =  90  -  axh 

The  above  equations  can  be  derived  with  the  use  of  Figure  1  below: 

NCP 


90-d 


h  =  LHA~    -  A.   -  a 


O 


.zi  =  uu-aih 


star 


Figure  1.    -  The  Astronomical  Triangle 


2.2   Lunar  Time 

As  we  have  defined  the  lunar  right  ascension  S3rstem  in  terms  of 
ascending  node,    it  would  only  be  appropriate,   that  we  define  the  lunar  time 
system  in  terms  of  the  draconic  month.      The  draconic  month  is  defined  by 
two  successive  passages  of  the  moon  through  the  same  node  of  its  orbit. 
The  draconic  month  is  equal  to  27.212220  mean  solar  days. 

The  expression  of  selenographic  longitude  in  terms  of  time  would  be 
as  follows: 


1  month 
1  day 
1  hour 
1  minute 
1  second 


27.212220 
.9070740 
.03779475 


mean  solar  days 
mean  solar  days 
mean  solar  days 


.0006299125     mean  solar  days 
.0000104985     mean  solar  days 


360 
12° 
30' 
30" 
0".5 


Figure  2  below  shows  the  relationship  in  draconic  times. 


local  mean 


astronomic  meridian 


eq.N 


LMDT  Lunar  mean  draconic  time 

LADT  Lunar  apparent  draconic  time 

RMDT  Reference  mean  draconic  time 

RADT  Reference  apparent  draconic  time 

eq.N  Equation  of  the  node 


Figure  2.    -  Draconic  Time  System 


From  Figure  2  it  is  evident  that: 

LMDT   +  eq.  N  =  LADT 
and  LMDT  -  RMDT  =    X 


where    X 


and    eq.N 


The  reduced  astronomic  longitude  with  the  positive 
sign  to  the  east 

The  lunar  equation  of  the  node  which  is  due  to  nutation 
(physical  libration  in  longitude  T  and  in  node  cr )  (see 
Appendix  C). 


From  the  actual  rotation  angles  of  the  moon,   we  find  that    eq.N  =T-o 

We  will  assume  in  this  section  that  the  meridian  is  not  affected  by 
polar  motion  and  that  the  rotation  of  the  moon  is  free  from  short  periodic 
irregularities. 

The  hour  angle  of  the  node  is  equal  to  the  lunar  right  ascension  of 
an  object  at  the  instant  of  its  upper  transit  over  the  meridian  from 
which  the  hour  angle  is  reckoned.      Therefore,   draconic  time  can  be 
determined  by  observing  transits  of  stars. 

The  lunar  time  intervals  are  almost  equal  to  the  terrestrial  time 
intervals  of  the  same  name  (91%  of  the  terrestrial  values),   which  is 
still  advantageous  from  the  psychological  point  of  view  because  a  change 
to  such  a  time  system  would  be  almost  unnoticeable  to  man.      The 
meridian  to  use  for  a  change  of  dates  could  be  the  180     meridian. 

2.3     Lunar   Pole  Stars 

As  stated  by  Gurevich    [5]  the  lunar  poles  describe  circles  in  the 
celestial  sphere  with  an  angular  diameter  of  3    .      This  circle  is  centered 
on  the  ecliptic  pole,   with  a  period  of  18.6  years.      Therefore,   first 
certain  stars,    then  others,   will  be  most  suited  for  use  as  pole  stars. 
Tables  of  altitudes  and  azimuths  must  be  compiled  for  the  stars  closest 
to  the  poles  at  a  particular  time.       The   mean   north   pole 
of  the  moon  described  in  the  earth  right  ascension  system  would  be 
approximately,     &  =   18h  ,    6    =   66°33T  .     We  see  from  the  FK4  catalogue 
that  the  following  stars  from  Table  1  are  near  the  north  pole. 
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Table  1.        Lunar  North  Pole  Stars 


Catalog  Number 

Magnitude 

659 

5.21 

701 

6.00 

664 

4.87 

685 

5.03 

a 


17h  32° 

68°  09' 

18h  36m 

65°  27' 

17'  37m 

68°46T 

18h  13* 

64°  23' 

and  approximately    ol  =   6    ,     6 
from  the  Boss  catalogue. 


-  66    33'  for  the  south  pole  in  Table  2 


Table  2.      Lunar  South  Pole  Stars 


Catalog  Number 

Magnitude 

a 

6 

6927 

5.3 

h         d 

5    32 

-  64°  15f 

7246 

4.5 

5h44° 

-  65°  45' 

7384 

5.2 

5h49B 

-  66°  55' 

7813 

5.8 

6h06B 

-  66°  02' 

7946 

4.9 

h          n 

6     11 

-  65° 35' 

With  a  knowledge  of  the  relationship  of  the  position  of  the  lunar  pole  to 
those  stars  the  astronauts  can  make  a  rough  north  alignment  with  a  probable 
accuracy  of  less  than  1  minute  while  using  only  a  gun  sight  type  sighting 
device. 
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2.4     Measurements  on  the  Lunar  Surface 

Presented  here  are  four  basic  types  of  measurements  which  can  be 
employed  on  the  lunar  surface  to  determine  astronomic  position. 

One  method  is  the  measurement  of  altitudes  or  zenith  angles  of  at  least 
two  stars  separated  in  azimuth,   combined  with  the  known  celestial  coordinates 
of  a  star  and  a  knowledge  of  time.      This  is  a  method  which  can  be  employed 
with  an  ordinary  sextant  or  theodolite. 

A  second  method  is  by  determination  of  a  relative  bearing  and  altitude  of 
a  star  combined  with  information  from  inertial  and  gravitational  sensors,   which 
is  the  method  used  by  the  Alignment  Optical  Telescope   (AOT)  on  the  lunar 
module. 

A  third  method  is  by  measurement  of  relative  azimuth  angles  between  at 
least  three  known  stars.      This  technique  can  be  used  with  an  ordinary  theodolite, 
but  also  lends  itself  well  to  the  use  of  an  automatic  instrument. 

A  fourth  method  is  by  star  photography,   which  can  use  either  the  first 
or  third,    or  other  methods  to  determine  position.      Here  we  can  use  a  photo 
theodolite,   or  other  metric  camera. 

In  the  following  section  these  four  methods  will  be  discussed  in  relation 
to  their  ease  of  employment  on  the  lunar  surface,    instrumental  adaptability 
and  accuracy,   and  other  pertinent  facts. 

2.4.1     Method  1:     Measurements  of  Altitudes 


The  measurement  of  altitudes  or  zenith  angles  in  order  to  determine 
astronomical  positions  can  be  accomplished  by  the  use  of  a  sextant  or 
theodolite.      The  obvious  advantage  of  this  method  lies  in  the  fact  that  both 
sextants  and  theodolites  are  comparatively  light  equipment,   which  could 
feasibly  be  included  on  a  lunar  mission.      The  sextant  is  a  familiar 
instrument  to  the  astronauts,   as  they  already  have  been  trained  in  its  use. 
Although,   the  training  required  to  allow  the  astronauts  to  effectively  use  a 
theodolite  could  be  accomplished  in  a  short  period  of  time.      The  big 
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advantage  in  using  this  equipment  is  that  astronomical  positions  can  be  taken 
at  several  points. 

As  stated  in  the  Geonautics  study  [7],   it  is  not  likely  though  that  a 
theodolite  would  be  practicable  for  the  initial  missions.      The  time  and  effort 
required  to  set  up  and  measure  is  far  out  of  proportion  to  the  many  other 
tasks  required  of  the  astronauts.      It  would  also  be  necessary  to  adapt  the 
fine  pointing  and  adjusting  screws  for  use  with  the  pressure  gloves  the 
operator  must  wear.      In  addition,   precise  sightings  through  the  helmet 
face  plate  would  be  most  difficult. 

A  hand  held  sextant  may  be  used  to  determine  a  position  on  the  lunar 
surface.      The  altitude  of  a  star  is  read  directly  from  the  scale.     With  the 
aid  of  sub-divisions  and  a  vernier,   altitudes  can  be  read  with  a  first-class 
instrument  to  one-tenth  of  a  minute  of  arc. 

The  problem  which  arises  with  the  use  of  a  simple  sextant  is  in  the 
determination  of  the  lunar  horizon.      In  order  to  eliminate  this  problem, 
we  can  consider  the  use  of  an  artificial  horizon  or  a  bubble  sextant,   so 
that  we  then  need  only  define  the  local  vertical  instead  of  the  horizon. 

The  artificial  horizon  is  a  rectangular  shallow  basin  of  mercury.      The 
observer  places  himself  so  that  he  can  see  the  image  of  the  body,   whose 
altitude  is  to  be  measured,   reflected  in  the  mercury.     He  points  the 
telescope  toward  the  artificial  horizon  so  that  he  uses  a  second  image  of 
the  heavenly  body  in  place  of  the  horizon.      By  moving  the  index  arm,   the 
two  images  are  made  to  coincide.      The  circle  reading  is  twice  the  altitude 
which  thereby  reduces  the  precision  of  the  instrument  to  12  seconds  of  arc. 

A  bubble  sextant  (bubble  octant)  [17]  contains  its  own  artificial  horizon. 
In  measuring  altitudes,   the  image  of  the  celestial  body  is  brought  into 
coincidence  with  the  image  of  a  circular  level  bubble  which  forms  the 
artificial  horizon.      The  altitude  is  then  read  on  a  micrometer  drum.      The 
fine  drum  is  divided  into  2  minute  intervals  which  can  be  interpolated  to 
the  nearest  1  minute.      In  a> bubble  sextant  such  as  the  Bendix  "Pioneer", 
provision  is  made  to  record  a  series  of  altitudes  automatically.     An 
average  reading  can  be  determined  from  these  readings. 
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These  instruments  are  not  capable  of  first  order  geodetic  accuracy. 
The  positional  error  can  be  a  maximum  of  50  meters  from  the 
sextant  with  the  artificial  horizon,   and  600  meters  for  the  bubble  sextant. 
This,   of  course,   does  not  take  into  account,    sextant  internal  errors, 
bubble  errors,   and  errors  in  the  star  tables.      In  addition,   the  recording 
of  each  measurement  using  the  common  sextant  will  burden  the  astronauts. 
With  both  instruments  each  star  must  be  identified  by  the  astronauts, 
leading  to  further  possibilities  of  error.      Time  recording  must  also  be 
done  manually. 

Using  this  method  (see  Figure  1)  we  can  find  azimuth  from  the 

formula  as  given  by  Mueller  [16] 

-  cos  d  sin  h 
sin  A  =  • 


sin  Z 

where    h  =  LHAq-X  -    a^ 

We  can  find  the  latitude  from  the  basic  formula  as  given  by  Mueller  [16] 
cos  Z  =  sind  sin^+cosd  cosh  cos$ 
and  we  can  determine  longitude  from  the  formula  as  given  by  Mueller  [16] 

cos  Z  -  sin  d  sin^* 


cos  h  = 


cos  d  cos^ 


2.4.2     Method  2:     Measurement  with  the  Alignment  Optical  Telescope 

A  relative  bearing  and  altitude  of  a  star  can  be  determined  by  the 
Alignment  Optical  Telescope  (AOT)  which  is  currently  installed  in  the  Lunar 
Module  (LM).     An  astronomic  position  on  the  lunar  surface  can  be  established 
by  the  use  of  this  instrument  in  conjunction  with  the  Inertial  Measurement 
Unit  (IMU)  and  the  Primary  Guidance  Navigation  System  (PGNCS)  accel- 
erometers.      The  great  advantage  of  this  method  lies  in  the  fact  that  no 
additional  equipment  or  training  of  the  astronauts  is  necessary  to  obtain  an 
astronomical  postion.     An  apparent  disadvantage  is  the  fact  that  only  one 
astronomical  position  can  be  determined  during  each  landing  due  to  the 
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fact  that  the  equipment  is  fixed  in  the  lunar  module  which  in  turn  is  fixed 
on  the  lunar  surface. 

The  AOT  as  described  in  the  Grumman  manual  [4]   is  a  unity  power 
periscope  type  device  with  a  60    conical  field  of  view.      Its  primary  pur- 
pose is  to  provide  a  method  of  establishing  an  inertial  reference  with 
respect  to  celestial  bodies  for  inflight  and  lunar  surface  operations.      It 
is  operated  manually  by  the  astronauts  and  is  mounted  on  the   PGNCS 
navigation  base  to  establish  a  mechanical  alignment  and  common  reference 
with  the  EMU.      The  AOT  has  a  moveable  shaft  axis  which  is  parallel  to 
the  LM  "X"  axis  and  a  Line  of  Sight  (LOS)  axis  which  is  the  center  of 
field  of  view  and  fixed  at  45    from  the  LM  +  "X"  axis. 

The  AOT  LOS  is  fixed  in  elevation  and  moveable  in  azimuth  to  six 
detent  positions.      These  detent  positions  are  located  at  60    and  are 
selected  manually  by  turning  a  detent  selector  knob  on  the  AOT.      Three 
of  the  positions  can  be  used  operationally  for  star  sightings  and  are  called 
forward,   left,   and  right  (see  Figure  3).      The  forward,   or  zero  detent 
position  places  the   LOS  in  the  x,z  plane  looking  forward  and  up.      The 
right,   or  +60    position  places  the  LOS  60    to  the  right  of  the  x,z  plane  and 
the  left  or  -60     position  places  the  LOS  60    to  the  left  of  the  x,z  plane. 
Each  of  these  positions  maintain  the  LOS  at  45    from  the  LM  +  X  axis. 
The  remaining  three  detent  positions  are  storage  positions  with  the 

o 

position  180    from  the  forward  position  being  the  normal  storage  position. 
The  reticle  control  makes  it  possible  for  manual  rotation  of  the  reticle 
for  use  in  lunar  surface  alignments.      The  AOT  reticle  pattern  consists  of 
cross  hairs  and  superimposed  archimedes  spiral  (see  Figure  4). 

The  vertical  cross  hair  is  an  orientation  line  designated  the  "Y"  line 
and  is  parallel  to  the  LM  X  axis  when  the  reticle  is  at  the  0°  reference 
position.      The  horizontal  cross  hair  is  an  auxiliary  line  designated  the  "X" 
line  and  is  perpendicular  to  the  orientation  line.     A  one  turn  spiral  is 
superimposed  from  the  center  of  the  field  of  view  to  the  top  of  the  vertical 
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FIGURE  3.    -  AOT   Line  of  Sight   Positions 
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cross  hair.      The  reticle  is  edge  illuminated. 

A  counter  located  on  the  left  side  of  the  AOT,    displays  an  angular 
readout  of  reticle  rotation.      The   counter  reads   in  degrees  to  within  ±0.02 
degrees  or  ±72  seconds.      Interpolation  is  possible  to  within  ±0.01  degrees. 

On  the  lunar  surface,   the  selected  star  can  be  sighted  in  either  the 
AOT  left, right  or  forward  positions.      The  astronaut  using  the  manual  drive 
knob,   adjusts  the  reticle  to  superimpose  the  orientation  line  on  the  target 
star.      This  reticle  angle  displayed  on  the  AOT  counter  is  then  manually 
inserted  into  the   Lunar  Module  Guidance  Computer   (LGC)  via  the  Display 
Console   (DSKY).      This  provides  the  computer  with  the  star  orientation 
angle   (shaft  angle).      The  astronaut  then  continues  to  rotate  the  reticle  until 
a  point  on  the  spiral  is  superimposed  on  the  target  star.      This  second 
angular  readout  is  then  entered  into  the   LGC  via  the  DSKY.      The  AOT 
detent  position  and  the  star  code  number  are  also  inserted  into  the  LGC 
via  the  DSKY.      The  LGC  can  now  calculate  the  angular  displacement  of  the 
star  from  the  center  of  the  field  of  view  b}r  computing  the  difference  between 
the  two  counter  readings.      Due  to  the  characteristics  of  the  reticle  spiral, 
this  A  angle  is  proportional  to  the  distance  of  the  star  from  the  center  of 
the  field  of  view.      Using  this  A  angle  and  the  proportionality  equation,   the 
computer  is  then  able  to  calculate  the  trunnion  angle.      During  the  lunar 
surface  phase  of  the  lunar  landing  mission,   the  Mission  Control    Center 
and  the  LGC  use  the  unit  gravity  vector  as  measured  by  the  LM  PGNCS 
accelerometers  to  align  the   LM  platform  and  to  estimate  the  latitude  and 
longitude  of  the  LM  on  the  lunar  surface. 

The  AOT  star  alignment  establishes  an  inertial  reference  by  employing 
the  following  formulas: 


SA  = 


(Si  X  S3)  X  Si 

Si  X  S2 


X 
Y 
Z 


MC  I 


[A]     [MCI] 
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Where  Si    and  S2    are  known  unit  vectors  of  stars  in  moon  centered 
inertial  space  and    MCI   is  the  unknown  coordinates  of  the  spacecraft  with 
respect  to  the  moon  centered  inertial  system. 


SB 


h 


Sim 
Sim  X  S2  m/X  SiM 

Sim  X  S2M 


U 


X 
Y 
Z 


SM 


[B]     [SM] 


Where  Sim    and    S2M    are  measured  unit  vectors  of  stars  in  stable 
member  space  of  the  stars,   and    SM    is  the  assumed  stable  member 
coordinates  referenced  to  the  moon  centered  inertial  coordinate  system. 

[Sa]     =     [Sb] 

[A]     [MCI]     =    [B]     [SM] 

[SM]     =     [B]"1    [A]     [MCI] 

['REFSMMAT]     =     [B]"1    [A] 
Therefore:  MCI     =     [REFSMMAT]1     [SM] 

Upon  comparison  with  the  gravity  vector  and  an  estimate  of  the  lunar 
radius  at  the  point,   a  position  on  the  lunar  surface  can  then  be  estimated. 

Considering  only  the  precision  of  the  counter  reader  (with  interpolation) 
we  see  that  the   positional  accuracy  may  be  as  poor  as   152  meters 
or  498  feet.      If  we  considered  the  errors  in  the  determination  of  the  unit 
gravity  vector,   and  the  errors  in  the  star  tables,   we  would  obviously  find 
that  this  system  is  not  adequate  for  determining  accurate  positions  on  the 
lunar  surface.      Considering  a  minimum  standard  error  of  20".  8,   we  will 
eliminate  this  system  from  consideration  for  geodetic  purposes. 
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2.4.3     Method  3:     Measurement  of  Relative  Azimuth  Angles 

Method  three  involves  the  simultaneous  determination  of  longitude, 
latitude,   and  azimuth  based  on  the  motion  of  stars  in  azimuth.     An 
ordinary  theodolite  can,    of  course  be  used,   but  it  again  suffers  the 
drawbacks  mentioned  in  Method  1. 

This  method  though,   lends  itself  to  a  concept  which  can  allow 
automatic  determination  of  astronomic  position  and  azimuth.     An 
"Electronic  Zenith  Camera"  (EZC)  is  currently  under  development  by 
the  Control  Data  Corporation  of  Minneapolis,    Minnesota.      It  is  described 
by  Carroll  as  a  portable  sensor-computer  combination  which  performs 
both  detection  of  stars  and  calculation  of  astronomic  position  [1]   [2]. 
Stars  are  detected  by  their  transits  across  a  fan  of  radial  slits  placed 
at  the  focal  surface  of  a  vertically  oriented  optical  system. 

The  principal  advantages  of  this  system  is  that  once  set  up,   the 
system  can  determine  its  astronomical  position  on  the  lunar  surface 
without  further  attention  from  the  astronauts.      The  system  itself  can 
be  easily  set  up  on  the  lunar  surface.      This  set  up  includes  approximate 
north  orientation,   and  approximate  levelling.      If  necessary,   the  system 
could  be  left  to  operate  after  the  astronauts  have  left  the  lunar  surface. 
Radio  transmission  of  current  updated  positions  could  be  relayed  to 
earth  if  necessary. 

It  is  considered  that  astronomical  longitude  and  latitude  can  be 
determined  to  within  0.3   (lcr)  seconds  and  azimuth  to  within  3   (la) 
seconds   (not  considering  star  catalogue  errors)  which  is  certainly 
acceptable  accuracy  at  this  time.      No  instrument  of  this  type  has 
currently  been  manufactured,  but  all  advance  planning  has  been  completed 
and  the  necessary  technical  skill  is  available  to  produce  this  instrument. 

The  EZC  can  be  broken  down  into  the  system  block   diagram  as  seen 
in  Figure  5. 
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Figure  5.    -  EZC  Block  Diagram 

We  will  examine  each  component  starting  with  the  sensor 
(see  Figure  6).  The    lens  system  is  an  f/3.3  system  having  a  3 

centimeter  aperture  and  a  10  centimeter  focal  length.      The  lens  system  has 
a  16  degree  total  field  of  view   (eight  degrees  maximum  off-axis  angle),    and 
can  form  an  image  of  a  point  source  restricted  to  a  tangential  dimension 
no  greater  than  10  seconds  of  arc  but  not  restricted  by  the  radial  dimension 
It  will  be  necessary  to  shield  the  lens  from  such  bright  bodies  as  the  sun 
and  earth. 

The  reticle  has  100  radial  slits   (at  the  focal  point  of  the  lens)  about 
1.4  centimeters  long  and  4.8  microns  wide. 
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Figure  6.    -  EZC  Sensor  Unit 
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Directly  behind  the  slit  reticle  is  placed  an  end-window  photomultiplier 
tube  which  collects  all  light  passing  through  the  slits.      The  output  of  this 
tube  is  amplified,   filtered,   and  processed  so  that  only  the  maximum 
signal  (peak  pulse)  will  be  recorded  and  timed.      The  time  input  can  be 
any  commercially  available  temperature  compensated  ciystal  oscillator. 

Looking  at  the  control  unit,    (see  Figure  7)  we  find  that  the  heart  of 
the  system  is  a  miniature  general  purpose  digital  computer  from  which 
the  sensor  is  controlled,   and  from  which  the  operator  receives  status  and 
calculated  information.     Also  included  are  a  precision  oscillator  (crystal 
oscillator),   a  power  source,   and  the  necessary  interface  electronics.      The 
level  reference  is  a  mercury  pool  autocollimator.      Miniature  electronic 
autocollimators  can  reflect  a  beam  off  the  surface  of  a  completely  enclosed 
mercury  pool  (dust  free  surface),   the  surface  of  which  is  covered  by  a 
silicone  fluid  which  damps  out  ripples.      If  the  unit  is  to  function  during 
the  lunar  night,    it  will  be  necessary  to  heat  the  mercury  pool,   due  to  the 
lunar  temperature  being  below  the  freezing  point  of  mercury. 

The  DC  reversal  motor  is  under  command  of  the  control  unit.      The 
motor  drives  the  instrument  about  the  vertical  axis  against  one  of  two 
stops  placed  180     apart.      The  reversal  bearing  is  composed  of  dissimilar 
metals  (brass  and  bronze)  thereby  eliminating  the  possibility  of  micro- 
scopic welds,   which  would  inhibit  accuracy. 

Prior  to  operation,   the  astronaut  must  level  the  instrument  to  within 
30  seconds.     Within  this  range  the  automatic  level  reference  supplies 
information  to  the  computer  so  as  to  measure  the  angle  between  the 
instrument  vertical  axis- and  the  true  local  vertical.      The  astronaut  must 
also  make  an  initial  north  (south)  -  alignment  of  the  instrument.      This 
can  be  accomplished  by  using  a  gunsight  mechanism  on  the  equipment  and 
sighting  with  the  aid  of  near  polar  stars  such  as  previously  described. 
This  wall  yield  an  initial  approximation  of  north,   with  the  exact  azimuth 
to  be  solved  for.  v 
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Figure  7.    -  EZC   Control  Unit 
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We  are  interested  in  detecting  stars  from  +  1.    to  +5.3  magnitude. 
It  has  been  determined  that  the  best  results  are  obtained  with  a  star 
transiting  at  least  one  slit  each  minute.      In  order  to  achieve  this  rate 
on  the  lunar  surface,   we  must  either  wait  approximately  27  times  longer 
for  a  solution  than  we  would  on  earth  (approximately  27  hours)  or  have  an 
internal  spinning  reticle  which  would  reduce  the  solution  time  but  would 
introduce  additional  errors  into  the  system  (approximately  1-arc  second). 

Star  identification  is  based  on  an  approximate  knowledge  of  the 
observer's  position,   and  north  alignment.      The  computer  determines  the 
expected  pattern  of  transit  times  from  the  stored  star  catalog.      This 
computed  pattern  and  the  measured  time  pattern  are  then  superimposed, 
and  these  target  pairs  which  correspond  closely  and  uniquely  in  magnitude 
constitute  tentative  identifications.      These  are  then  used  to  adjust  the 
assumed  position  and  azimuth  and  thereby  to  adjust  slightly  the  computer 
time  pattern  so  as  to  obtain  a  better  fit  with  the  measured  pattern.      From 
here  the  final  solution  sequence  preceeds  by  solving  a  linearized  version 
of  the  following  equation  derived  from  equations  2.2: 

Q  ft 

cot(  A0  +  /3j  )  sin(LHA~  -  X  -  a   i  )  -  sin$cos  ( LHA~  -  X  -  a  t  ) 

+  tan  di  cos  $  =  0 


where 


A0  -       The  azimuth  of  the  Oth  slit  (dependent  on  north  alignment) 

jSj  -       The  fixed  angle  between  the  Jth  slit  and  the  Oth 

slit 
LHA^   -       Apparent  lunar  reference  meridian  hour  angle  of  the  node 

(LADT) 
X  -       Longitude 

0{  -       Lunar  right  ascension  of  the  i  th  star 

<P  -       Latitude 

di  -       Lunar  declination  of  the  i  th  star 
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This  linearized  version  is  then 


0  =  sin  Aj  cot  Z  d?+  cos*(  tan*-  cosAj  cot  Z  )  dh  -  dAc 


where  Z  is  the  computed  zenith  distance 


ft 


h  =  LHA~  -X-ttj 

Aj  =  A0  +  P< 


and 


cot  Z  = 


(  sin  d  sin*  +  cos  d  cos  h  cos^)  sin  Ai 
-  cos  d  sin  h 


The  measured  quantities  are  azimuth  and  time.      Latitude  and  longitude 
of  the  station  are  assumed. 


Since  there  are  more  measurements  than  unknowns,   we  may  use  a 
least  square  solution,  with  a  Kalman  filtering  technique  [10]  based  on  the 
Newton-Raphson  iteration  formula,    so  that  matrix  inversions  will  not  be 
required.      This  iterative  procedure  is  based  on  the  following  equation: 

Ym      -  y     -  r»    nf      [  aN  XN  +  LN  ] 

Xn+i  -An    -cNaN    r  -    ,  rt     _      t, 

[  1  +  aw  CN  aN  ] 


where: 


matrix  Xn    =    3X1 


~d*~ 

m       » 

dh 

x0   = 

0 

_dAo_ 

0 

I—        - 

matrix  aN    =    ia3   =    (  sin  Aj  cot  Zi    -    cos*  (tan*- cos  Aj  cot  Zi)  -  1) 

ft 

Ln    =    cot(A0  +  Aj)  sin(LHA^  -   X  -  G^  )  -sin*  cos(LHAq 

-  X-  a.  t  )     +    tan  di  cos* 
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and  matrix  CN+i  =  CN  - 


[  1  +  aN  CN  aft  ] 
where  matrix       C0     =  3C3  =(AI  £)~ 
where  the  rows  of  A  are  the  rows  of  ai  to  as  . 

We  can  see  that  the  incorporation  of  a  new  measurement  does  not 
require  a  new  matrix  inversion.      Only  one  inversion  is  required  and  that 
is  for  Co  which  is  only  an  inversion  of  a  3  x  3  matrix.     After  the  first 
five  marks  no  old  observation  information  need  be  saved.      The  values  of 
the  elements  of  matrix  C  are  indicative  of  the  quality  of  the  results. 
Matrix  C  is  the  covariance  matrix,  but  only  under  the  assumption  that 
there  are  no  systematic  errors  in  the  residual  vector  L  .      There  are 
systematic  errors  in  the  residual  vector  L  ,   due  to  the  timing  error 
(  T  libration  )  and  the  errors  in  star  right  ascension  and  declinations. 
In  order  to  cancel  alignment  errors  in  the  optical  -  reticle  system  and 
normality  errors  between  the  instrument  axis,   the  instrument  is  rotated 
180   ,   and  the  astronomic  position  is  again  determined.      The  average  of 
these  positions  is  then  used.      The  level  status  of  the  reversal  bearing 
is  also  measured  and  it  is  applied  as  a  final  correction. 

2.4.4     Method  4:     Measurement  with   Photography 

The  determination  of  astronomical  position  by  photography  can  be 
accomplished  with  a  photo-theodolite  or  other  metric  camera    using  methods 
such  as  one  or  three.      The  drawbacks  of  a  photo- theodolite  are  the  same 
as  for  an  ordinary  theodolite  and  thus  would  tend  to  make  it  impractical 
for  early  missions. 

As  described  in  the  Geonautics  study  [7]  ,   a  metric  camera  could  be 
used  as  a  zenith  camera  or  otherwise  fixed  in  elevation  to  photograph  stars. 

26 


Various  sophisticated  systems  such  as  gravity-operated  automatic  leveling 
devices  with  gimballed  camera  support,   electronically  operated  automatic 
leveling  devices,  gyroscopic  devices,  etc.,   could  be  considered  to  level  the 
camera.      The  simplest  and  most  feasible  device  would  seem  to  be  the 
conventional  spirit  bubbles.      These  bubbles  could  be  placed  on  a  camera 
support  such  as  a  tripod,    so  that  the  camera,   when  placed  on  the  support, 
could  be  theoretically  leveled  within  the  range  of  two  orthogonally  placed 
precise  level  vials  inside  the  camera  body.      These  precise  levels  can  be 
photographed  during  each  exposure,   thus  providing  a  permanent  record  of 
the  level  reference.      Drawbacks  in  this  system  tend  to  lie  in  the  fact  that 
the  bubbles  are  affected  by  differential  heating  from  external  light,   so  that 
they  must  be  protected  both  from  radiant  heat  and  from  light.      The  bubble, 
of  the  sensitivity  required  in  the  camera,   would  be  termed  a  "slow"  bubble 
in  a  conventional  surveying  instrument.      It  would  be  even  slower  on  the 
lunar  surface  due  to  the  moon's  smaller  gravitational  force. 

This  method  provides  a  permanent  and  true  record  of  field  data  as 
film  can  be  removed  from  camera  and  returned  to  earth.      Most  advantageous 
is  the  fact  that  a  great  number  of  measurements  can  be  taken  in  a  minimum 
amount  of  time  after  the  initial  set  up.      There  is  a  requirement  for  extra 
equipment  such  as  a  tripod,   but  it  could  be  retractable  and  made  of 
aluminum  or  some  similar  material       Timing  of  each  photo  would  also  be 
required,   but  this  could  be  done  by  an  internal  clock,   which  would  furnish 
the  time  to  be  included  in  the  edge  data. 

2.5     Error  Analysis 

Let  us  first  consider  the  errors  which  contribute  to  the  lunar  apparent 
seleno equatorial  coordinates.      In  the  transformation  of  the  mean  geoequatorial 
coordinates  to  the  mean  ecliptic  coordinates,   the  relatively  small  errors  in 
ol  ,    6  ,    and  £  can  be  considered  negligible.      Therefore,  /3  and  ^  can  be 
considered  errorless,      hi  the  transformation  of  the  mean  ecliptic   coordinates 
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to  the  mean  selenoequatorial  coordinates  we  consider  that  the  error  in  CI 
is  negligible,  but  there  is  definitely  an  error  in  inclination  (I).      Let  us 
make  a  conservative  estimate  of  20"  for  this  error.      Differentiating 
equations   (2. 1)  we  see  that  the  results  are  similar  to  that  given  by 
Mueller  [16]  . 

CI 
ddm  =  cos  a  a    sin  I  d  (  X  -  CI )   +  (  sec  da  sec  £  cos  I  -  tan  dm  tanp* ) 

CI 
dp  +  sin  a  m  dl 

but  d  ( X  -  CI )  =  0 

dj3  =  0 

CI 
therefore     ddE  =  sina  0  dl  and  therefore  the  maximum  value  of 

ddB  =  sin  a  B   d  I     =     20".  0  sina   B 

a  2 

and       da,    =  (sec  dB  cos  I  -  secda  tanda  sin/3)  d  (X  -  CI ) 

ci       0  .  a 

-  cos  a  b   sec  p  sec  da  sin  I  d  p  -  tan  d0  cos  a  a    d  I 

but  again     d  (X  -  CI )    -   0 
dj3  =   0 

CI  Ci  CI 

therefore     da  m     =    -tandn  cosa  B    dl    =   -20".  0  tandm  cos  B 


Considering  the  transformation  of  mean  into  apparent  selenoequatorial 
coordinates,   we  will  use  the  following  equations: 

,  „,  Ci  CI         d  Ai  ai  n       3  Ai  ai    ,         3  Ai  ai 

da^pp  =  da  B     +  -r^ — dcr+    r-j        dl  +    .  j    x  dda 

OCT  dl  ddn 

.        3Aiai         Q     SBibi  3Bibi  3  Bx  bi    _, 

+      da-    da»  +     dp    *  +  TTdI  +Td7~dda 


m 


dBibi        CI      dCiCi   .     CI      SCiCi     1T         SCiCi    „ 

+  ¥a^d^+-^5da»  +~dT—  dI  +^d7~ddB 

3Didi         0      dDx  di     ,T  \  dDi  di     5J 

+  -d^¥dam+~it"^dI   +-TZ~dd'» 


(2.3) 
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and  a  similar  expression  for     ddsapp,    considering  in  both  cases  a 
negligible  error  in  proper  motion. 

In  the  following  differential  expressions  we  will  assume: 

cos  I  =1 

dm  =     Limited  to  45 

p  =     1' 

a  =    i« 


The  differential  quantities  will  be  assigned  the  following  conservative 
values: 


da 

= 

10" 

= 

.000050  radians 

dp 

= 

10" 

= 

.000050  radians 

dl 

= 

20" 

= 

.000097  radians 

ddn 

= 

20" 

= 

.000097  radians 

dT 

da  „ 

= 

15" 

20" 

= 

.000014  radians 
,000097  radians 

} 


(estimate  from  [13] ) 


(estimate) 


,  a, 

Aiai  =  (P0t  +  a  )  (cos  I  -  sin  I  sin  a  D  tan  dB  ) 


d  Ai  ai 
Sct 


da   =  (cos  I  -  sin  I  sin  a:  0  tan  dB  )  da 
=   (1  -  .02676  sina  0   tandm)da    =    da 
=  10". 0 


5  Ai  ai 
b  I 


d  A 

dI=(P0t+a)(  cos  I  sina  B    tan  dm  -  sin  I )  dl 

=  (-.0015059)  (-.97324)dl 

=  0".0 


5Ai  ai 

3  dB 


d  ft  2 

ddm    =    (  P0 1  +  a)  (sin I  sina  B     sec   dm)  ddffl 

=   (-.0015059)  (.05352)  ddB 

=  0".0 
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Al^  da  ,       =       (Pot  +a)  (-sin  I  cosaa  tand=)  dot  m 
3a  u  CI 

B  =       (-.0015059)    (- .02676)  da  B 

0".0 


CI 
Bibi  =  p  cos  a  B  tan  dm 


3Bi  bi  fi 

dp     =  cos  a  B  tan dm  dp 


Bp 


10".  0  cos  a?  tan dffl 


^^dl       =  0".0 

dl 


3Bi  bi  ^2 

— — —  ddm  =  pcosa  m  sec   dm  ddm 
o  dm 

=  .00058  dda 
=  0".0 


SBibi   ,    CI  ft         ,     ,    CI 
7T  da  B        =      psina  „  tandn  da  B 

Sam  Ct 

-  (.00029)  dam 

-       0".0 


o  fl 

Ci  ci  =  K  sin  (Lap   -  U  +  180  )cosI  cos  a  csecdmdl 


— ^r^dl    =-K  sin  (Lap   -  ft  +  180°)  sinl  cosa  B    secdmdl 
o  I 

=  (.000097)(1)(.02676)(1)  (1.4)  dl 
=  0".0 

—  ddm    =  K  sin  (LqP  -  ft  +  180  )  cos  I  tand=  secdD  cosa  „  ddB 

O  dm 

=   (.000097)  (1)    (.02676)   (1)   (1)    (1.4)  dl 
.     =    0".0 
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ir-  da  a        =         K  sin  (Lap  -  &+  180  )  C  =  SI  sin  a  B  secdB  da 

=        (.00097)    (1)    (.02676)    (1)    (1.4)  da. 


0".0 


o  O 

Di  di  =  K  cos  (Lap   -  £i    +  180  )  sin  a  B    secdB 


3  Di  di 

— ~ — di     -  0".0 

~T~r~  ddm   =  K  cos  (Lap   -  ft  +  180°)  sin  a  B  tandm  secdB  ddm 
a  dB 

=  (.000097)   (1)   (1)   (1)   (1.4)  ddB 
=  0".0 

BDidi,    ^        =       Ksin  (  LaP  -  ft  +  180°)  cos  a  „    secdmdaB 

— ry   Q<-J-    m  "i 

da  *  =        (.000097)  (1)  da« 

0M.0 


We  wall  now  consider  the  standard  error  in  lunar  right  ascension.     We 
will  use  the  quadratic  form  of  equation  2.3,    eliminating  terms  which  are 
zero.     We  will  assume  that    da*«  and 

dBibi    . 

— > dp 

dp      "** 

are  dependent  terms  while  all  other  non-zero  terms  are  independent. 

Therefore, 

(d«if      =       (da?)2  +  (^-di)3    +  (^-dpf 

+  ^da?dp 
3p 

=  400  tan2dmcos2a  B  +  100 
+  100  tan2dmcos2a9 
-  200  tan2dDcos2a9 


Q,  [~  2  2     Q 

and  da  sap?  =   V 300  tan  dm  cos  a  B   +100 

31 


A       f  a  ft 

Ai  ai  =  +  (  P0  t  +  cr )  sin  I  cosa  B 


BAi  ai    .  .    _         „  ft    , 

— r^ —    dcr    =  sin  I  cosa  _   da 

ft, 

=  .02676  cosaBdCT 

Max.  value  =  0".5  (will  be  neglected) 


dAi  ai  ,    *  ft 

— — - —  d  I    =  ( p0t  +  a  )  cos  I  cos  a  m  d  I 
a  I 

=  -.0015059  dl 

=  0M.0 


^Ai  ai  , , 
r*      ddB    =  0M.0 

O  affl 


5Aiaxf   ,    ft         #tAjl  .    ,    .    _.ft,..ft 


-75—  da  a  =  -(Pot  +cr)  -sin  I  sina_  da  m 


(-.0015059)  (-.02676)   (1)    (da  ■) 


Bi  bi  =  psina  m 


-d  Bi  bi  d  p=  _   sina  n   dp 
dp 

^idl     =0".0 
o  I 


BBibl 


3t-^ddm    =  0".0 


BBib,1  ,     ft  ft,    ft 

— q    da,       =  pcosamda„ 

bav  0 

=  -.00029  da  V 

=  0".0 


Ci  c[  =  -K  sin  (Lap  -  ft +  180°)  cos  I  sina  B   sindm 

-K  sin  (Lap-  ft  +  180°)  sin  I  cos  dm 
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— Jr^dl     =  (  +  K  sin  (Lap  -  Cl  +  180°)  sin  I  sin  a  B   sin  dm 

ol 

-  K  sin  ( Lap  -  Cl  +  ISO    )  cos  I  cos  da )  d  I 
=  -.00007  dl 
=  0".0 


— r — ~  ddm   =  ( -  K  sin  ( LaP  -  Cl    +  180°)  cos  I  sin  a  a  cos  dm 

odm 

+  K  sin  (Lap  -  Cl  +  180°)  sin  I  sindm)  ddB 
=  -.00007  dda 
=  0".0 


Xq   d«  B         =       -Ksin(LaP  -  Ch+  180°)  cosIcosaffl  sindB 

B  n 

.00007  day 

=       0".0 


Di  di  =  -  K  cos  (  LaP  -  Cl  +  180°)  cos  a  a    sindB 

^^dl     =  0».0 

->  \    1  ddm   =  -  K  cos  ( La  P  -  Cl  +  180  )  cos  a  n  cos  dn 

=     .00007  ddB 
=  0".0 


3DidLT   ,    Cl  _  0  Q 

7f   dQ„       -       Kcos  (Lap  -  Cl+  180  )  sina  B  sindB 


9aV 


.00007  da*? 
0".0 
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We  will  now  consider  the  standard  error  in  lunar  declination.     We  will 
use  the  quadratic  form  of  the  similar  declination  equation  2.3,    eliminating 
terms  which  are  zero.     We  will  assume    ddn    and 


dBibidp 
dp 


are  dependent  terms. 


Therefore, 


(ddsappf  =  (ddBf  +  r|j^dp)a 

+     ^  dpddB 

dp 


=    400  sin2  ol  n    +    100  sin3  ol  m 

3       & 

-200  sin  ol  m 

2       ft 

=    300  sin  ol  m 


a, 

and  ddsapp      =     17".3sinaB 


Thus,   we  can  see  that  if  stars  with  near  zero  declination  are  used,   we 
can  effectively  limit  the  error  in  lunar  right  ascension.      Let  us  examine 
the  circumstances  if  we  use  stars  with    ±10    declination  and  ±  10    right 
ascension.      The  maximum  error  will  then  be: 

dasapp    =V300  (.03109)  (1)  +100  =  10"5 
ddsapp      =    17". 3  (.17633)     -    3".l 

Thus,   we  can  see  that  there  can  be  an  optimum  time  for  which  to  obtain 
a  position.      If  we  consider  limiting    0£"B"    to  values  of  350    right  to  010 
and  170    right  to  190  ,   we  will  immediately  limit  the  error  in  declination 
to  a  maximum  value  of  3".l.      See  Appendix  E  for  stars  which  fall  into 
this  area. 
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Next  we  will  consider  errors  in  time.     We  have  defined 

LMDT   +  eq.N  -  LADT 

or 

RMDT  +  eq.N  =  RADT 

We  can  consider  that  RMDT  is  errorless  so  that  the  error  then  occurs 
only  in  eq.N  which  is  dependent  upon    cr   and    T.      So  that  we  may  say 

dT2  =  da2  +  dT2  +  dadr 


.$ 


or  dT  =  10  . 0    draconic  time  or  5".0  of  arc. 

We  will  now  consider  the  standard  error  of  each  method  in  deter- 
mination of  astronomic  position. 

For  Method  1,   we  will  first  present  the  differential  formula  for 
latitude  [16]  : 

d  3?    =     -  sec  A  d  Z    -  cos  §  tan  A  d  h 

-  sin^cosec  Z  sec  d  -  tan  d  cot  Z  dd 
cos  A 

We  desire  the  azimuth  to  be  near  0  or  180  degrees  so  that 

dS*    =     -  dZ  -  (sin^cosec  Z  sec  d  -  tan  d  cot  Z  )  dd 

and  if    Z    equals  approximately  45  ,   then  for  a  station  at  or  near  the  equator 
d    will  also  equal  about  45      so  that 

d<£    =     -  dZ  -  dd 

Therefore,   for  Method  1,   these  errors  are  independent,   letting    dZ  =  7".0, 
and  considering 

(d§)2     =     (7".0)2  +  (17". 3  sinaa  f 

and  d$         =     7".0  to  18".  7 

Thus,   an  error  in    h    will  have  little  effect  on  latitude  determination.      The 
error  in  latitude  determination  is  basically  related  to  the  approximate  right 
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ascension.      The  differential  for  azimuth  [16]  is 

dA  =  (sin^cosec2Z  -  sin  d  cosec  Z  cot  Z  )  dh 
+  cos <i>sin  h  cosec2Z  dd  +  cot  Z  sin  Ad$ 

If  we  consider  the  zenith  distance  near  90    we  see 

dA  =  (sin^>-  sind)  dh  +  cos<£sinhdd 

If  the  station  is  near  the  equator,     <£  will  be  close  to  0  .      Therefore, 

dA  =  -  sin  ddh  +  sin  hdd 

Finally,   we  will  select  a  pair  of  stars  close  to  the  meridian,   preferably 
one  with  90    north  declination  and  one  90    south.      Therefore,     dA  =  0 
for  stations  near  the  equator  when  we  combine  the  pair.      If  we  use  only 
a  sextant,   we  will  have  no  means  of  accurately  turning  the  angle  between 
a  star  and  landmark  so  that  no  permanent  azimuth  could  be  computed. 
The  equation  for  longitude  [16]  is 


where 


a. 


a 


sapp 


h 
T 

ATn 

Therefore    dX 


X    =   ttsapp    +   h    -    (T    +    ATn) 


the  apparent  right  ascension 
the  hour  angle  as  determined 
the  local  time  on  time  keeper 
the  chronometer  correction 


da 


sapp 


+  dh  -  dTn  -  dAT 


and 


dh  =     -  (cosec  A  sec<£dZ   +  cot  A  sec ^d^) 

-  (tan^cosec  Z  sec  d  -   tan  d  cot  Z  sec  d)  dd 
sin  A 


The  error  in    h    may  be  eliminated  by  observing  at  an  azimuth  of  90    or 
270    and  by  observing  east  and  west  stars  in  pairs  which  are  symmetrical 
to  the  meridian  and  at  the  same  altitude.      Considering    tan<i>  is  zero,   the 
star  pairs  will  eliminate  the  declination  errors  and  zenith  error  terms. 

36 


So  that  dX  =  da,™  -  dT-dAT, 


sapp 


d  ATn    has  little  effect  on  the  moon  for  any  ordinary  chronometer,    so  that 

d*=  da^  -  dT 


Therefore  dX2  =  (daSappf  +  (dT)    -  dTdaBPP 

2       ^ 

-  300  tan2  dffl  cos  a.  B  +  100  +  25.0  - 
5.  V300  tan2  dffl  cos2  a  n    +  100 
Limiting    dB    to  45°,   we  see 

dX  =  8". 6  to  18". 0 

For  Method  2,   we  have  already  stated  that  the  low  precision  of  the 
AOT  precludes  further  consideration. 

For  Method  3,   with  impersonal  operation  and  high  equipment  precision, 
we  can  consider  that  the  only  errors  presented  will  be  systematic  in  nature, 
and  these  specifically  will  be  the  systematic  errors  in  star  coordinates 
and  time.      Thus,   we  can  consider  that  an  astronomic  position  by  the  use 
of  this  method  will  yield  errors  near  the  equator  (declination  <  10  ) 

d<±>=  ddgapp  =  3".l  (for  optimum  time) 
dX  =  dajjpp  -  dT 

dX2=  (da6app)2+(dT)2  -  dT   daJPp 
dX   =  9".l 

and  dA  =  sin  Ai  cot  Zd$+  cos<3?  (tan<£>  -  cos  Aj  cotZ)dX 

near  the  equator,   we  say 

d  A  =  (sin  Aj  d<±>  -  cos  Aj  dX)  cot  Z 

The  largest  Z  we  can  hope  for  is  8°,   and  assigning  an  azimuth  of  the  jth 
slit  to  be    45     ,   we  find: 
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Therefore, 


dA    =  25  d<£2  +  25  dX2    (considering  no  dependence) 
dA    >   50". 0 


For  Method  4,   the  errors  are  basically  the  same  as  for  Method  1, 
if  a  panoramic  camera  is  used,    and  basically  the  same  as  for  Method  3, 
if  the  zenith  camera  is  used  (considering  camera  operating  at  optimum 
time).      See  Table  3  for  comparison  of  results. 


Method 
1 
2 
3 

4* 


d$ 


7".0  -  18". 7 


>    20". 8 


dX 


S".6  -  18M.0 


>    20". 8 


3".  1  (optimum  time)       9".l 
3".l  -  18".  7  8"- 6  -  18".  0 


dA 

not  directly 
obtainable 
not  directly 
obtainable 

>  50". 0 

>  0".0 


*Dependent  upon  horizon  or  zenith  camera  type,    not  including 
instrument  or  observation  errors. 


Table  3.    -  Estimated  Standard  Errors  (Station  Near  the  Equator) 


2.6     Conclusions 

Due  to  the  possibility  of  limiting  the  error  in  star  coordinates,    it 
would  seem  that  an  automatic  instrument  as  described  in  Method  3,   would 
be  the  one  which  would  be  able  to  take  best  advantage  of  this  point. 
Missions,   of  course  cannot  be  scheduled  so  that  these  star  coordinates 
could  be  utilized.      It  would  be  only  coincidence  if  the  current  short 
duration  landings  occured  at  the  proper  time.      Thus,   the  advantages  of 
an  instrument  which  could  determine  position  automatically  after  the 
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astronauts  depart  is  apparent.      During  future  missions,    if  the  landings  are 
near  the  lunar  equator,    it  would  seem  that  the  Electronic  Zenith  Camera/ 
Zenith  photography  would  take  best  advantage  of  the  stars  with  near  zero 
declination.      If  landing  were  far  removed  from  the  equator,   optimum 
times  would  need  to  be  recomputed  due  to  the  increased  error  in. lunar 
right  ascension,   which  would  result  from  the  higher  star  declinations. 
The  format  of  the  current  missions,   therefore,   dictates  that  we  use  an 
automatic  instrument  which  could  take  advantage  of  stars  which  offer 
minimum  coordinate  errors.      The  ease  of  set  up  and  minimum  training 
required  for  the  astronauts,   combined  with  the  inherent  accuracy  of  the 
system,    show  that  the  Electronic  Zenith  Camera  or  similar  instrument  is 
the  best  possible  choice  of  equipment  for  the  current  missions.     An 
accurate  azimuth  can  only  be  determined  by  the  use  of  an  horizon  camera 
at  this  time. 

It  must  again  be  stated  that  the  Electronic  Zenith  Camera  has  not  yet 
been  manufactured,   and  could  not  be  included  in  immediate  missions. 
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3.   DEFLECTION  OF  THE  VERTICAL 

Deflections  of  the  vertical  relating  the  normals  of  the  astronomic 
and  selenographic  systems  can  be  obtained  by: 

(1)  Direct  astronomical  measurement  of  latitude  and  longitude 
with  subsequent  comparison  with  the  selenodetic  coordinates 
extended  from  a  datum  by  triangulation  or  traverse.      The 
results  are  affected  by  accumulation  of  error  in  the 
triangulation  or  traverse. 

(2)  Integrating  gravity  anomalies  over  large  areas  surrounding 
the  point  where  the  deflection  is  required. 

(3)  Computing  the  deflections  of  the  vertical  due  to  topography, 
and  then  refining  this  system  using  some  model  of  isostatic 
compensation. 

(4)  Computing  the  deflections  of  the  vertical  from  known  spherical 
harmonics. 

The  first  method  depends  on  an  assigned  deflection  at  a  datum,  thus 
is  of  no  use  in  this  paper.      The  second  method  requires  gravity  anomalies 
over  large  areas,   and  these  of  course  are  not  available  except  by  obtaining 
them  through  method  four.      But,   method  four  is  of  little  value  in  deflection 
computations  as  the  harmonics  reflect  moon-wide  rather  than  local 
distribution.      The  computation  of  deflections  of  the  vertical  requires  more 
local  data  near  the  station.      Thus,   we  will  investigate  method  three,   but 
will  use  method  four  as  a  check  on  the  magnitude  of  the  deflections. 

In  addition,   the  effect  of  other  bodies  on  the  deflection  of  the  vertical 
will  be  calculated. 

3.1     Computation  of  Topographic  Deflections  of  the  Vertical 

Following  the  method  of  Jordan-Eggert  [11]  adapted  for  the  moon,   we 
consider  the  influence  of  the  visible  masses  in  the  neighborhood  of  a  station. 
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For  this,   we  imagine  the  terrain  around  the  station  divided  into  compart- 
ments by  means  of  concentric  circles  and  radial  lines,    so  that  we  obtain 
a  great  number  of  vertical  columns  between  a  given  reference  surface 
and  the  moon's  physical  surface. 

If  the  components  of  this  deflection  in  the  meridian  and  the  prime 
vertical  are  expressed  by  A£    and  A77   respectively,   then  we  have 


3  0  H  P  " 

A£    =    — (sina2    -    sinax)   loge   r3  /  n 

4  fe  IT  ma 

-  3  0  Ho" 
kr\     -       .  r.  -  K      (cosa2   -    cos  a  1)  loge  r2/ri 
4  c'  it  m  a 

where       0  a  -  The  surface  layer  density 

S  -  The  mean  density  of  the  moon 

H  -  Height  above  reference  surface 

a  -  The  radius  of  the  moon 

aia2  -  The  compartment  azimuths 

ri  ,  r2  -  The  outer  and  inner  radius  of  compartments 

The  reference  surface  which  will  be  used  is  defined  as  a  moon  centered 
sphere  of  radius  1,737.988  kilometers. 

We  will  use  the  following  values  for  the  constants  as  given  in  [6]. 

8  =  2.670  gm/cm3     (assumed)* 

0  ffl  =  3.343  gm/cm3 

a  =  1,737,988  meters 

p  =  206,265" 


*  This  figure  was  obtained  from  averaging  figures  in  various  sources 
for  a  surface  layer  of  10KM.  An  error  in  this  figure  will  produce 
a  proportional  error  in  the  deflections. 
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and  we  obtain 


A£      =     0".0226  H(sino:3   -    sinax)  loge  r2/r! 
A77      =     0".0226  H(cos0£3  -    cosa  x)  loge  r2  /  r x 
Note:     H    in  meters 


These  two  values  are  to  be  computed  for  each  terrain  section  around 
the  point  0,   and  then  the  two  components  of  the  deflection  of  the  vertical 
for  the  point  0  are  equal  to  the  algebraic  sum  of  all  values  of  A  £  or  of 
A 77  as: 

m  m 

i    =     S    A£j  77    =    S     A^i 

1=1  1=1 

We  have  so  far  assumed  that  the  whole  of  the  attracting  stratum  is 
in  the  horizon  of  the  station,   which  for  the  moon  would  be  grossly  in 
error  when  only  a  short  distance  from  the  station.      Therefore,    it  is 
necessary  in  order  to  secure  the  required  degree  of  accuracy  to  take 
into  account  the  curvature  of  the  moon's  surface.      The  effect  of  this 
curvature  is  to  throw  the  compartment  surface  below  the  horizon  of  the 
station.      Therefore,   when  significant,   we  will  use  the  equation  developed 
by  A.    R.    Clarke  ,  given  in    [11]  . 


A  4  =  0".0226  H  (sina3  -  sinax  )  log, 


tan  1/4  0  a 


tan  1/4  0  1 
+  cos  l/2  02   =  cos  1/2  0  1 


A  77  =  -  0".0226  H  (cosa  2  -  cosa  1)  log 


tan  1/4  0  1 
+  cos  l/2  02  -  cos  1/2  0  1 

in  which  62  is  r2  reduced  to  an  arc 
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-.  r2  ( in  km  ) 

©2     -    .   „»„  „^       radians 
1,737.988 

and    ©i    is    ri    reduced  to  an  arc 

©i     =     * ^   radians 

1  1,737.988 

Any  values  for  either  the  constant  r2  /  ri    or  for  the  constant  difference 
(sinG:2   -    sin^i)  may  be  adopted  arbitrarily.     Advantage  is  taken  of  this 
fact  to  adopt  such  values  as  would  make  the  computation  as  simple  and 
rapid  as  possible.      The  values    r2  /  n  =  1.19362053  and  (sinQ!2     -    sinQ:i)  = 
0.25  were  adopted.     With  these  values,   the  deflection  equations  will  be: 

A£         =     .001"  h  (in  meters) 

A77         =     .001"h  (in  meters)     (the  reference  line  is  shifted  90  ) 

As  a  result  then,   of  this  particular  selection  of  arbitrary  constants 
defining  the  limits  of  the  compartments,   the  deflection  produced  at  the 
station  by  the  material  lying  above  the  reference  surface  in  any  compart- 
ment is  expressed  as  one  second  of  arc  for  each  kilometer     the  mean 
elevation  is  above  the  reference  surface. 

The  arbitrary  selection  of  the  adopted  values  of  (sino;2  -  sinO!].) 
and  r2  ,  ri    was  guided  by  three  considerations.      First,    it  was  important 
to  be  as  rapid  as  possible  in  determining  the  contribution  of  each  section. 
The  compartments  must  be  small  enough  to  bring  the  accuracy  of  the 
method  within  required  limits  and  yet  large  enough  to  avoid  excessive 
amounts  of  detail  in  the  computation.      The  compartments  should  be 
compact  areas,   not  long  and  narrow,    in  order  to  facilitate  the  estimation 
of  the  mean  elevation  within  each  compartment. 

The  radii  of  the  adopted  circles  separating  the  compartments  are 
shown  in  Table  4.  Each  group  of  sixteen  compartments  (four  in  each 
quadrant)  forms  a  ring.      The  rings,    for  convenience  of  designation, 
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were  assigned  serial  numbers  commencing  with  the  inner  ring,    (see 
Figure  8 ) . 

It  was  determined  that  the  curvature  of  the  moon  had  to  be  taken 
into  account  at  ring  46.      The  outer  radii  of  all  succeeding  rings  were 
computed  iteratively  on  the  IBM  360-75  computer,    (see  Appendix  A). 

It  was  found  by  Hayford  [8]  ,   that  on  the  earth,    the  computed 
topographic  deflections  are  much  larger  than  the  actual  deflections  of  the 
vertical,    in  that  an  influence  is  in  operation  which  produces  an  incomplete 
counterbalancing  of  the  deflections  produced  by  topography,   leaving  much 
smaller  deflections  in  the  same  direction.      The  deflections  of  the  vertical 
are  due  to  irregularities  in  the  distribution  of  the  masses  composing  the 
earth.      These  irregularities  may  occur  either  as  a  result  of  irregularities 
in  the  surface  of  the  earth   (topography),    or  as  a  result  of  irregularities  in 
the  distribution  of  the  densities  beneath  the  surface. 

The  distribution  of  the  densities  below  the  surface  of  the  earth  is 
invisible  and  largely  unknown,   but  the  theory  of  Isostasy  postulates  pre- 
cisely such  a  relation  between  sub-surface  densities  and  surface  elevations. 

If  the  earth  were  composed  of  homogeneous  material,    its  figure  of 
equilibrium,   under  the  influence  of  gravity  and  its  own  rotation,   would  be 
an  ellipsoid  of  revolution. 

The  earth  though  is  composed  of  heterogeneous  material  which  varies 
considerably  in  density.      If  this  heterogeneous  material  were  so  arranged 
that  its  density  at  an}'  point  simply  depended  upon  the  depth  of  that  point 
below  the  surface,    or,   more  accurately,    if  all  the  material  lying  at  each 
equipotential  surface  (rotation  considered)  was  of  one  density,   a  state  of 
equilibrium  would  exist  and  there  would  be  no  tendency  toward  a  rearrange- 
ment of  masses. 

If  the  heterogeneous  material  composing  the  earth  were  not  arranged 
in  this  manner  at  the  outset,   the  stresses  produced  by  gravity  would  tend 
to  bring  about  such  an  arrangement,   but  because  the  material  is  not  a 
perfect  fluid,   the  rearrangement  will  be  imperfect.      In  the  partial 
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RING 


Outer  Radius   (KM) 


RING 


Outer  Radius   (KM) 


1 

.0080 

36 

3.9211 

2 

.0095 

37 

4.6803 

3 

.0114 

38 

5.5865 

4 

.0136 

39 

6.6682 

5 

.0162 

40 

7.9592 

6 

.0194 

41 

9.5003 

7 

.0231 

42 

11.3398 

8 

.0276 

43 

13.5354 

9 

.0330 

44 

16.1562 

10 

.0393 

45 

19.2844 

11 

.0470 

46 

23.0184 

12 

.0561 

47 

27.4754 

13 

.0669 

48 

32.7956 

14 

.0799 

49 

39.1461 

15 

.0953 

50 

46.7267 

16 

.1138 

51 

55.7757 

17 

.1358 

52 

66.5781 

18 

.1621 

53 

79.4741 

19 

.1935 

54 

94.8707 

20 

.2310 

55 

113.2545 

21 

.2757 

56 

135.2082 

22 

.3291 

57 

161.4304 

23 

.3928 

58 

192.7602 

24 

.4688 

59 

230.2080 

25 

.5596 

60 

274.9951 

26 

.6680 

61 

328.6057 

27 

.7973 

62 

392.8567 

28 

.9517 

63 

469.9958 

29 

1.1359 

64 

562.8462 

30 

1.3558 

65 

675.0283 

31 

1.6184 

66 

811.3203 

32 

1.9317 

67 

978.2927 

33 

2.3057 

68 

1185.5146 

34 

2.7522 

69 

1448.1099 

35 

3.2850 

70 

1793.1006 

Table  4.    -  Ring  Outer  Radius 
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Reference 


FIGURE  8.    -     Lunar  Template 
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rearrangement  some  stresses  will  still  remain  as  different  positions  of 
the  same  horizontal  stratum  may  have  somewhat  different  densities,   and 
the  actual  surface  of  the  earth  will  depart  from  the  ellipsoid  of  revolution 
in  the  sense  that  above  each  region  of  deficient  density  there  will  be  a 
bulge  or  bump,   and  above  each  region  of  excessive  density  there  will  be  a 
hollow.      The  bumps  on  the  earth  would  be  mountains  and  the  hollows 
would  be  ocean  floors.      The  excess  material  represented  by  that  portion 
of  the  mountains  above  the  reference  surface  will  be  compensated  for  by 
a  defect  of  density  in  the  underlying  material  and  vice-versa  for  the 
oceans.      This,   of  course,    can  be  explained  by  the  famous   Pratt-Hayford 
theory  of  Isostasy. 

If  we  now  consider  the  moon,   we  see  that  the  consequences  of 
hydrostatic  theory  as  stated  by  Kopal  [15]  demand  that  given  a  sufficiently 
long  time  the  material  which  on  a  short  time-scale  may  behave  as  a 
solid  is  bound  to  get  crushed  under  its  own  weight  to  settle  to  a  form  of 
minimum  potential  energy,   which  is  a  sphere.      The  departures  of  the 
lunar  globe  from  spherical  form  permissible  for  a  figure  of  equilibrium 
in  the  prevailing  field  of  force  are  indeed  small.      To  this  extent,   the 
moon  may  seem  to  be  in  hydrostatic  equilibrium.      However,   this  agree- 
ment becomes  noticeably  worse  when  we  compare  the  momenta  of  the 
lunar  globe  with  their  theoretical  values.      Finally,  by  considering  the 
mechanical  ellipticity,   we  see  that  its  value  is  so  much  at  variance  with 
the  requirements  of  hydrostatic  equilibrium  that  the  discrepancy  must  be 
considered  as  real. 

Therefore,    in  spite  of  the  approximate  prevalence  of  hydrostatic 
equilibrium  in  the  lunar  interior,   as  attested  by  the  nearly  spherical  form 
of  its  globe,   the  moon's  motion  around  its  center  of  gravity  reveals  the 
presence  of  small  but  unmistakable  departures  from  this  equilibrium. 
On  the  other  hand,    O'Keefe  [19]  was  able  to  show  by  comparison  between 
calculated  Bouguer  anomalies  and  the  essentially  free  air  anomalies  of 
Muller  and  Sjogren,   that  the  topography  of  the  moon  could  be  isostatically 
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compensated  and  that  to  a  rough  first  approximation  the  moon  is  in 
isostatic  equilibrium.      He  considers  an  actual  crustal  layer  on  the 
moon,    for  which  he  states  the   Pratt-Hayford  theory  of  Isostasy  could 
be  applied. 

Under  the  assumption  of  Pratt's  hypothesis  of  Isostasy,   we  will 
introduce  a  surface  of  compensation  at  the  depth    T    below  the  reference 
surface.      If  we  denote  the  changes  of    A£    and    A  77   due  to  the  changes  in 
topography  by    A£0    and  At?0  ,   then  the  sums    A£    +   A£0    and    A17  +  A770 
are  the  terms  which,   when  summed,   give  the  total  deflection  of  the 
vertical.      From   Jordan- Eggert  [11]  : 


A?o    .     -3&eTp"  .    SInOl)l0g,  \* 

4  77  0B  A 


3A0   TP"     .         „.  „.   .,         /ra  +  V  T2  +  r22 


At70  = — (C0S0!2    -  cos Q!i) log 


e 


vri  +  V  T  +  ri 


A  0f    =     8H 


Therefore: 


A£°  +  A£    =    ;  -7J P  "      (sina2         sinai)    [loge-^ 

4  7T    C*  a     A 


ri 


V  T3  -f  r23 


r  +  W  +  n2 
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Then  for  a  definite  value  of  T  we  can  set  up  a  table,  from  which  for 
every  value  rx  (because  r2  =  constant  Xrx)  the  factor  F  can  be 
taken.      Then  we  will  have 

A£    +    A£0    -    fA4 
At)   +    At70    =    fAtj 

O'Keefe  states  [19]  that  for  20%  difference  in  lunar  crust  and  lunar 
mantle  density,   a  depth  of  compensation  of  25  kilometers  is  in  order, 
while  for  10%  difference,   a  depth  of  compensation  of  50  kilometers  is 
in  order.      Further  extending  these  assumptions,    calculations  for  depths 
of  compensation  of  25,   50,   75  and  100  kilometers  has  been  made  using 
the  program  and  subprogram  in  Appendix  A. 

Table  5  gives  the  various     F    factors  for  each  ring  owing  to  a 
specific  assumed  depth  of  compensation. 

Calculations  have  also  been  made  for  125,   250,   500,    1,000  and  1,500 
kilometers  depth  of  compensation  primarily  to  obtain  a  full  set  of  values. 

There  is  direct  evidence  of  excess  mass  in  the  maria  over  what 
would  be  called  for  by  Isostasy.      The  explanation  of  this  excess  mass 
which  was  discovered  by  Muller  and  Sjogren  of  JPL  is  subject  to  various 
interpretations.      O'Keefe  [19]   states  that  this  excess  mass  may  be  due 
to  the  maria  being  filled  with  some  sort  of  material  which  arrived  after 
the  formation  of  the  highland  topography.      H.    C.    Urey  [21]  states  that 
this  mas  con  effect  could  be  produced  by  a  flat  circular  slab  of  chrondritic 
material.      For  instance,    Urey  calculated  that  in  Mare  Imbrium,   this 
disc  would  be  some  4  kilometers  thick  and  be  670  kilometers  in  diameter 
(mass  =  5  X  10      grams). 

J.    G.    Stipe  [20]  states  that  these  mas  cons  could  be  caused  by  meteorites 
which  have  bombarded  the  moon  and  have  been  of  necessary  size  and  at 
necessary  depth  to  create  these  mass  anomalies. 
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RING 


25  km 


50  km 


75  km 


100  km 


2 

0.9997 

0.9998 

0.9999 

0.9999 

3 

0.9996 

0.9998 

0.9999 

0.9999 

4 

0.9995 

0.9998 

0.9998 

0.9999 

5 

0.9994 

0.9997 

0.9998 

0.9999 

6 

0.9993 

0.9996 

0.9998 

0.9998 

7 

0.9992 

0.9996 

0.9997 

0.9998 

8 

0.9990 

0.9995 

0.9997 

0.9998 

9 

0.9988 

0.9994 

0.9996 

0.9997 

10 

0.9986 

0.9993 

0.9995 

0.9996 

11 

0.9983 

0.9991 

0.9994 

0.9996 

12 

0.9979 

0.9990 

0.9993 

0.9995 

13 

0.9975 

0.9988 

0.9992 

0.9994 

14 

0.9971 

0.9985 

0.9990 

0.9993 

15 

0.9965 

0.9983 

0.9988 

0.9991 

16 

0.9958 

0.9979 

0.9986 

0.9990 

17 

0.9950 

0.9975 

0.9983 

0.9988 

18 

0.9941 

0.9970 

0.9980 

0.9985 

19 

0.9929 

0.9965 

0.9976 

0.9982 

20 

0.9915 

0.9958 

0.9972 

0.9979 

21 

0.9899 

0.9949 

0.9966 

0.9975 

22 

0.9879 

0.9940 

0.9960 

0.9970 

23 

0.9856 

0.9928 

0.9952 

0.9964 

24 

0.9828 

0.9914 

0.9943 

0.9957 

25 

0.9795 

0.9897 

0.9932 

0.9949 

26 

0.9755 

0.9878 

0.9918 

0.9939 

27 

0.9708 

0.9854 

0.9903 

0.9927 

28 

0.9651 

0.9826 

0.9884 

0.9913 

29 

0.9584 

0.9792 

0.9861 

0.9896 

30 

0.9504 

0.9752 

0.9834 

0.9876 

31 

0.9408 

0.9704 

0.9802 

0.9852 

32 

0.9294 

0.9646 

0.9764 

0.9823 

33 

0.9158 

0.9578 

0.9718 

0.9789 

34 

0.8996 

0.9496 

0.9664 

0.9748 

35 

0.8804 

0.9399 

0.9599 

0.9699 

Table  5.    -     F    Factors 
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Table  5.    -     F    Factors   (Continued) 


RING 


25  km 


50  km 


75  km 


100  km 


36 

0.8577 

0.9283 

0.9521 

0.9641 

37 

0.8309 

0.9145 

0.9429 

0.9571 

38 

0.7994 

0.8981 

0.9319 

0.9489 

39 

0.7627 

0.8787 

0.9188 

0.9390 

40 

0.7200 

0.8556 

0.9032 

0.9272 

41 

0.6712 

0.8285 

0.8847 

0.9133 

42 

0.6163 

0.7965 

0.8627 

0.8966 

43 

0.5558 

0.7593 

0.8368 

0.8769 

44 

0.4908 

0.7161 

0.8063 

0.8535 

45 

0.4232 

0.6669 

0.7707 

0.8260 

46 

0.3557 

0.6115 

0.7293 

0.7936 

47 

0.2910 

0.5505 

0.6818 

0.7559 

48 

0.2319 

0.4852 

0.6282 

0.7122 

49 

0.1803 

0.4175 

0.5687 

0.6624 

50 

0.1370 

0.3501 

0.5045 

0.6065 

51 

0.1022 

0.2858 

0.4373 

0.5451 

52 

0.0751 

0.2272 

0.3695 

0.4795 

53 

0.0545 

0.1763 

0.3040 

0.4118 

54 

0.0392 

0.1338 

0.2435 

0.3445 

55 

0.0280 

0.0996 

0.1902 

0.2806 

56 

0.0199 

0.0731 

0.1452 

0.2225 

57 

0.0141 

0.0530 

0.1086 

0.1722 

58 

0.0099 

0.0381 

0.0800 

0.1304 

59 

0.0070 

0.0272 

0.0582 

0.0969 

60 

0.0049 

0.0193 

0.0418 

0.0710 

61 

0.0035 

0.0136 

0.0299 

0.0513 

62 

0.0024 

0.0096 

0.0212 

0.0368 

63 

0.0017 

0.0067 

0.0150 

0.0262 

64 

0.0012 

0.0047 

0.0105 

0.0185 

65 

0.0008 

0.0033 

0.0074 

0.0130 

66 

0.0006 

0.0023 

0.0051 

0.0091 

67 

0.0004 

0.0016 

0.0035 

0.0063 

68 

0.0003 

0.0011 

0.0024 

0.0043 

69 

0.0002 

0.0007 

0.0016 

0.0029 

70 

0.0001 

0.0005 

0.0011 

0.0019 
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Wise  and  Yates  [22]   state  that  mascons  are  produced  by  mantle  plugs 
upwelling  into  giant  impact  bases  punched  through  the  lunar  crust,  followed 
by  volcanic  filling  of  the  remainder  of  the  crater  above  the  plug. 

No  matter  which  theory  is  accepted,   the  direct  result  will  be  a  further 
deflection  of  the  vertical.      The  area  where  the  mascon  is  located,   will 
deflect  the  vertical,   proportionately  larger  than  the  height  of  this  area 
would  indicate.      The  exact  numerical  increase  in  the  deflection  would  be 
dependent  upon  the  mass  distribution  as  offered  in  the  theories  above. 
We  could  consider  these  mass  increases  as  additional  point  masses  and 
compute  their  effect  on  the  deflections  accordingly,   or  we  could  select  a 
station  which  would  minimize  their  effect. 

Let  us  consider  a  central  station  which  is  quite  distant  from  the  major 
mass  anomalies  so  that  their  effect  will  be  insignificant  (very  low  F  factor) 
in  the  overall  result.      For  instance,    considering  any  station  near  the 
equator  from  20    east  longitude  westward  to  40    west  longitude  will  limit 
the  contribution  of  these  mass  anomalies  on  the  deflection  of  the  vertical 
as  the  station  will  be  greater  than  600  kilometers  from  the  outer  edge  of 
these  anomalies.     We  can  see  for  a  station  greater  than  600  kilometers 
from  a  large  mass  anomaly  that  the  largest  F  factor  is  only  about  .  015 
(down  to  100  km    compensation  depth)  and  the  contribution  of  the  mass 
anomalies  are  then  generally  negligible  on  the  overall  result.      For 
topographic  deflections  only,   the  area  of  these  mass  anomalies  will 
generally  fall  into  one  compartment  only,   and  will  in  turn  affect  the 
result  very  little.      For  this  reason,    I  have  selected  the  center  of  the 
crater  Landsberg  "A"  as  a  starting  point.      Its  coordinates  are  328    57'W 
and  0     13' N.      This  point  was  also  selected  due  to  availability  of  charts. 
Due  to  the  nearly  concentric  area  in  terms  of  height  of  the  center  of 
Landsberg  "A"  and  the  use  of  a  1:250,000  chart,    I  have  neglected  rings 
1  to  35,   which  are  approximately  the  inner  3  kilometers  of  rings.      Rings 
36  to  50  were  computed  using  the  AMS  series  topographic  lunar  map  sheet 
1  edition  1  of  3  -  64.      This  map  is  a  modified  stereographic  projection 
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with  a  point  of  origin  4,410.10  kilometers  from  the  point  of  tangency  of 
the  map  plan  and  a  selenographic  sphere  with  a  radius  of  1,737.988 
kilometers.      The  contour  interval  is  250  meters  based  on  a  vertical 
datum  centered  at  Mosting  "A"  at  an  elevation  of  7,000  meters.      The 
horizontal  datum  was  also  based  on  the  center  of  Mosting  "A"  which  was 
considered  to  be  latitude  3°  10'  47"  S  and  longitude  354°  50'  13".      Rings 
51   to  63  were  computed  using  an  edition  1  -  AMS  sheet  2  topographic 
lunar  map,    scale  1:2,000,000  with  the  same  datum,    etc.,   as  above. 
Rings  64  to  68  were  computed  using  an  edition  2  -  AMS  sheet  2 
topographic  lunar  map,    scale  1:5,000,000  with  the  same  datum,   etc., 
as  above. 

Rings  beyond  68  were  neglected  due  to  their  approach  to  the  lunar 
limb;  and  their  coinciding  reduction  in  accuracy  of  height  determination. 

Tables  6  and  7  give  the  heights  of  each  compartment  for  computations 
in  the  prime  vertical  and  in  the  meridian  respectively. 

Table  8  shows  the  computed  deflections  for  assumed  depths  of 
compensation. 
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Compensation  Depth 

4 

25 

km 

+ 

2". 934 

50 

km 

+ 

3". 332 

75 

km 

+ 

2". 234 

100 

km 

+ 

0".883 

125 

km 

- 

0".547 

250 

km 

- 

6". 118 

500 

km 

- 

16". 792 

1,000 

km 

- 

30". 184 

1,500 

km 

- 

38". 529 

CD 

L'04".450 

n 

+  10". 983 
+      20". 709 

+  27".  757 
+  33". 123 
+  37". 447 
+  49". 548 
+  1' 11". 984 
+  11  38".  947 
+  1' 47". 427 
+  2f  17".480 


Table  8.    -  Deflections  Considering  Depths  of  Compensation 


We  cannot  consider  that  the  value  of    ©    will  remain  constant  down  to 
the  depths  which  we  are  considering.      Thus,   as  depths  increase,   the 
average  density  for  this  surface  layer  will  vary.     We  consider  that  the 
average  density  for  these  depths  are  as  follows  in  Table  9   (interpolated 
from  Kaula  [12] ) . 
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KM  Average  Density  F 1 

25  2.670  1 

50  3.000  1.1236 

75  3.140  1.1760 

100  3.210  1.2022 

125  3.250  1.2171 

150  3.263  1.2221 

250  3.336  1.2494 

500  3.343  1.2521 

1,000  3.300  1.2360 

1,500  3.343  1.2521 


Table  9.    -  Average  Surface  Density  for  Selected  Depth 


The  F  1  factor  in  Table  9  is  a  deflection  multiplier,   and  its  effect 
on  deflections  is  given  in  Table  10. 
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Compensation  Depth 


r\ 


25  km 

50  km 

75  km 

100  km 

125  km 

250  km 

500  km 

1,000  km 

1,500  km 


+ 

2' 

.934 

+ 

4' 

'.118 

+ 

2' 

.627 

+ 

1' 

.062 

- 

0' 

.666 

- 

7' 

'.644 

- 

21' 

.025 

- 

37' 

.307 

- 

48' 

'.242 

— 

1'04' 

.450 

+  10 
+  23 
+  32 
+  39 
+  45 
+  1'01 
+  1'30 
+  2'02 
+  2' 14 
+  2'17 


".983 
".269 
".642 
".820 
".447 
".905 
".131 
".298 
".509 
".480 


Table  10.    -  Deflections  Considering  Changes  in  Density  for 
Changes  in  Depth  of  Compensation 


The  graph  in  Figure  10  shows  the  relationship  of  depth  of  compensation 
to  deflections,   considering  no  change  in  surface  density,   and  considering  a 
change  in  the  surface  density  for  changes  in  the  depth  of  compensation. 
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The  maximum  topographic  deflections  on  the  moon  will  most  likely 
occur  at  a  highland-mar ia  dividing  line.      For  instance,   a  north-south 
deflection  was  examined  at  the  southern  wall  of  the  crater  Pitatus  at 
latitude  31°  10'  S  and  longitude  346    10'.     With  a  few  assumptions   (due 
to  lack  of  large  scale  charts  on  the  area)  it  was  determined  that  this 
deflection  was  approximately  5  minutes.      It  is  quite  possible  that  other 
larger* deflections  exist,   but  probably  will  not  greatly  exceed  this  value. 
It  is  interesting  to  note  that  the  deflections  may  be  greatly  minimized  by 
choosing  a  station  which  is  surrounded  by  nearly  concentric  rings  of  equal 
elevation.      This  will  yield  nearly  zero  deflections  for  the  inner  zones. 
This  ideal  station  could  be  at  the  center  of  a  concentric  crater,   at  the 
top  of  an  isolated  rise  on  flat  ground,    or  at  a  central  point  of  a  flat 
area.      If,   of  course,   a  station  were  selected  which  did  not  satisfy  this 
concentric  condition,   the  inner  rings  must  be  considered. 

3.2       Computation  of  Spherical  Harmonic  Deflections  of  the  Vertical 

In  order  to  compare  the  values  obtained  by  the  method  of  isostatic- 
topographic  deflections,   we  will  use  the  formula  as  given  by  Dragg  [3]  to 
compute  the  deflections  by  the  use  of  spherical  harmonics. 


00 

-  1  N  r  "NT    4-    1 

i  = S      £      (AnmcosMX  +  BNMsinMX)     — —  PNM  (sincp) 

yCOSCp      f^g    fpo  LN   -    1 

(N  -  M  +  1)  „  ,    .       1 

sin<p  -    —      _  1 "■  PN    +i  ,  ■   (sincp)j 


1  ~        1        M       .  N  x 

*7  =  : ■     ^       *->    — — r  (ANMsmM  X  -  B^cosMA)  pNM  (sirup) 

ycos  cp     N^  2    ^o     N  - 1  v  '         v  ' 


Dragg  [3]  neglected  the  20,    21,   and  40  spherical  harmonic  terms  in 
his  computation  for  the  earth,  but  they  have  been  included  in  this  paper. 
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I  will  use  the  spherical  harmonics  coefficients  given  by  Blackshear  [6], 
(see  Appendix  D).      These  are  CNm  and  Smm  coefficients  which  must  be  con- 
verted to  ANM  and  BNm  above. 

By  definition  we  see  that     CNm   =   Jnm   and    SNM   =   Knm. 

For  Heiskanen  and  Moritz   [9],  we  see 

A'NM  =  GM  rN  JNM 
B'NM    =  "GM  rN  K^ 


and  that  from  Dragg  [3] 


_  (N  -  1) 

Anm    -  ,  A  NM 

rN   +2 

_    (N  -   1) 

.dnm   -  a  mm 

N   +S 

r 


Thus  we  see  that  in  terms  of  coefficients     Cm    and  Snm 


CD 

1  N 

i  =  S     S     (CnMcosMX    +    SNMSinMX)    |(N  +   l)PNM(sin<p) 


COS<P        F2K4 

sin(p  -  (N  -  M  +  1)  PN  +  ]#  M  (sin<p)] 


00 

-1 


rj  -  S      2      M(CNMSinMX    -    SjmcosMX)   PNM(sin<p) 

COSCP         m=2    m=o 


This  equation  was  evaluated  for  the  Landsberg  "A"  station  using  the 
IBM  360-75  computer,   the  program  of  which  is  contained  in  Appendix  B. 
The  results  of  this  evaluation  are 

£    =    -    12'.  184 
77    =     2'13".707 

These  values  fall  between  the  isostatic-topographic  results  of  250  km 
depth  of  compensation,   and  the  results  of  topographic  -compensation  only, 
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when  surface  density  changes  are  not  considered.      These  values  fall 
between  the  isostatic-topographic  results  of  250  kin  and  1,500  km  depths 
of  compensation,   when  surface  density  changes  are  considered. 

3.3     Effect  of  Other  Bodies  on  the  Deflection  of  the  Vertical 

It  is  obvious  that  the  attractive  effect  of  the  sun  and  the  earth  will 
affect  the  deflection  of  the  vertical  on  the  lunar  surface.      The  following 
developments  refer  directly  to  the  earth  but  apply  equally  as  well  to  the 
sun.      Planets  also  have  an  effect  but  for  all  practical  purposes  this  effect 
can  be  considered  insignificant. 

Following  the  method  of  Jordan  -  Eggert  [11]  adapted  for  the  moon,   we 
see  that 


Deflection     = 


3  Me  A    sin  Z  cos  Z 

MB  r 


where 


and  where 


Me 

Mm 

A 

Z 

r 

Me 

Mffl 


-  The  mass  of  the  earth 

-  The  mass  of  the  moon 

-  The  lunar  radius 

-  Zenith  Distance  of  the  earth 

-  Earth-moon  distance 

=  81.301 


and 


A° 

~" 3" 

r 


(1737. 63 f 
(384,400) 


9.27  X  10 


-8 


Therefore: 

Deflection 


2.253  X  10     pM  sin  Z  cos  Z 
4". 646912274  sin  Z  cos  Z 
2". 323  sin  2  Z 
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For  the  effect  of  the  sun  we  will  use  the  formula 


3  Ma  A    sin  Z  cos  Z 

Deflection     =    — — — 5 

Mn  r 


Where  ~-    =     27,158,000 


M 


B 


A3  (      1,737.63       f  is 


Deflection     =     0".026  sin  Z  cos  Z 
=     0".013sin2Z 

which  is  negligible  in  comparison  with  other  uncertainties  introduced. 

We  will  now  divide  the  deflection  of  the  vertical  caused  by  the  earth's 
attraction  into  a  northward  and  westward  directed  component.      If    A    is 
the  azimuth  of  the  earth  counted  from  the  north,   then  we  have  for  the  two 
components  of  the  deflection  of  the  vertical 

£e     =     2".  32  sin  2  Z  cos  A 

?7e     =     2".  32  sin  2  Z  sin  A 

These  deflections  are  not  constant  for  a  station  but  vary  in  accordance 
with  the  change  in  earth-moon  distance  and  the  geometric  librations  in 
latitude  and  longitude.      They  also  vary  with  the  physical  librations  but 
these  are  insignificant. 

3.4     Conclusions 


Due  to  lack  of  knowledge  of  the  lunar  interior  and  its  isostatic 
properties,   we  cannot,   with  any  degree  of  confidence,    state  the  magnitude 
of  the  deflections  of  the  vertical  at  a  given  point. 

For  the  computation  of  the  topographic  deflections,   we  used  lunar 
charts  with  a  lunar  centered  spherical  surface  as  the  equipotential 
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reference  surface,   and  from  these  the  compartment  relative  heights  were 
determined.      If  we  consider  the  effect  of  an  offset  sphere  on  this  computation, 
we  can  see  immediately  that  the  relative  heights  would  change  in  accordance 
with  the  translation  of  the  sphere's  center.      For  instance,    if  the  center 
of  the  reference  sphere  were  offset  toward  the  earth,   the  relative  heights 
of  the  eastern  compartments  calculated  in  this  paper  would  be  consider- 
ably less,    and  the    r\   deflections,   considerably  smaller.      Thus,   we  will 
need  better  height  determination  on  the  moon ,   before  the  topographic 
deflections  can  be  considered  to  be  sufficiently  accurate.      Due  to  the 
relative  ease  of  computation  of  the  spherical  harmonic  deflections,    it  may 
be  best  to  use  them  at  present  instead  of  the  topographic  deflections, 
especially   in  areas  which  have  local  concentric  heights.      The  computed 
deflections  could  be,   of  course,   corrected  for  the  effect  of  the  earth. 

If  in  the  future,    Method  (1)  is  employed,   and  relative  deflections 
between  stations  were  obtained,   we  would  be  in  a  much  better  position  to 
predict  the  absolute  deflections  at  these  same  stations  using  the  topo- 
graphic method. 

The  employment  of  Method  (2)  would,    if  extensive  enough,  yield  the 
absolute  deflections  of  each  station.      It  is  considered  that  in  this  method, 
lies  the  best  possible  results.      But,   alternately,   the  most  time  and  effort 
would  be  required  over  the  other  methods. 

Further  work  can  be  accomplished  using  Method  (3).      For  instance, 
Karl  Jung  in  Report  No.   41  of  The  Ohio  State  University,    Institute  of 
Geodesy,    Photogrammetry  and  Cartography,   has  given  the  formulas  to 
obtain  approximate  values  of  the  influence  of  lunar  craters  on  the 
deflection.      These  computations  could  become  particularly  useful  for 
stations  offset  from  the  crater  centers,   as  they  would  yield  a  rapid 
evaluation  of  the  deflection  attributed  to  the  inner   rings  of  the  template. 
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APPENDIX  A 
F  FACTOR  COMPUTATION  PROGRAM  AND  ASSOCIATED 
MOON  CURVATURE  SUB-  PROGRAM 

The  main  program  computes  the  "F"  factor  for  a  given  ring,  given 
the  depth  of  compensation.    The  associated  sub-program  computes  the 
change  in  outer  ring  radius  due  to  moon  curvature.    The  main  program  list 
of  variables  is: 

IR  Indexing  variable  for  ring  number 

J  Indexing  variable  and  integer  compensation  depth 

Al  Inner  radius  of  ring 

A2  Outer  radius  of  ring 

T  Real  number  compensation  depth 

FA  Compensation  depth  squared  plus  outer  radius  squared 

FB  Compensation  depth  squared  plus  inner  radius  squared 

FC  Outer  radius  plus  the  square  root  of  FA 

FD  Inner  radius  plus  the  square  root  of  FB 

FE  FC  divided  by  FD 

FF  The  natural  log  of  FE 

FG  The  outer  radius  divided  by  the  inner  radius 

FH  The  natural  log  of  FG 

F  The  F  factor  equal  to  1  -  FF/FH 

The  sub-program  Mooncu  uses  an  iterative  procedure  employing  the 
secant  method  to  compute  the  outer  ring  radius  taking  moon  curvature  into 
account. 

The  sub-program  list  of  variables  is: 

THETA1  The  inner  radius  divided  by  the  lunar  radius 

expressed  in  radians 

THETA2  The  outer  radius  uncorrected  for  lunar  curvature 

expressed  in  radians 

THETB2  An  approximate  outer  radius  corrected  for 

curvature  expressed  in  radians 
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Tl 

T2 

T2B 

T21B 

T21 

Til 

T3 

T3B 

SOC 


ASOC 


THETC2 


THETA1/4 

THETA2/4 

THETB2/4 

Tangent  of  T2B 

Tangent  of  T2 

Tangent  of  Tl 

T21/T11 

T21B/T11 

The  natural  log  of  T3  minus  the  cosine  of  § 

THETA1  plus  the  cosine  of  |THETA2  minus 

the  natural  log  of  the  outer  radius  divided  by 

the  inner  radius  which  is  given  by  definition. 

This  is  considered  the  first  approximation. 

The  natural  log  of  T3B  minus  the  cosine  of  \ 

THETA1  plus  the  cosine  of  §THETB2  minus  the 

natural  log  of  the  outer  radius  divided  by  the 

inner  radius  which  is  given  by  definition. 

This  is  considered  the  second  approximation. 

By  the  secant  method  where: 


X 


i+l 


=  Xfcj  f(X,)  -  X,  i(Xi-,) 
f(Xi)  -f(X1_1) 


csoc 

BSOC 


we  see 

rmn-rrp-  -  THETA2(ASOC)  -  THETB2  (SOC) 
THETC2  -  ASQC  _   SQC 

The  difference  between  X1+1  and  X^ 
Outer  radius  in  kilometers 
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APPENDIX   B 

LEGENDRE   POLYNOMIAL   COMPUTATION   PROGRAM  AND 
SPHERICAL  HARMONIC   DEFLECTION   PROGRAM 


The  program  computes  the  Legendre  polynomials  and  the  deflections 
of  the  vertical  in  both  the  meridian  and  the  prime  vertical  using  the  13th 
degree  porynomial  coefficients  as  given  b}'  Blackshear  [6]  for  the  position 
of  the  center  of  Landsberg  "A".      The  Legendre  polynomials  were 
computed  to  the  fourteenth  degree  as  required  and  the  subscript  15 
represents  the  Oth  element  [P(15,15)  =  P(0,0)].      The  program  list  of 
variables  is: 


P(15,15) 
C(14,14) 

S(14,14) 

ALAT 

ALONG 

U 

T 

I,J,KA,M 

MA,MB,N 

PSI 

APSI 
PNU 

APNU 


An  array  which  contains  the  Legendre  polynomials 

An  array  which  contains  the  Blackshear  C 

coefficients 

An  array  which  contains  the  Blackshear  S 

coefficients 

Latitude  of  computation  point 

Longitude  of  computation  point 

Sine  of  latitude 

Cosine  of  latitude 

Indexing  variables 

Indexing  variables 

Evaluation  of  a  single  term  of  the  deflection 

in  the  meridian 

Sum  of  PSI  Terms 

Evaluation  of  a  single  term  of  the  deflection  in 

the  prime  vertical 

Sum  of  PNU  terms 
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PPSI  Multiplication  of  APSI  by  constants  for  final 

value  of  deflection  in  meridian 

PPNU  Multiplication  of  APNU  by  constants  for  final 

value  of  deflection  in  prime  vertical 
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APPENDIX  C 
PHYSICAL  LIBRATION  EQUATIONS 

As  stated  by  Kolaczek  [14] ,   physical  libration  is  the  physical 
deviation  of  the  moon's  motion  from  that  described  by  Cassini's  Laws. 
It  can  be  treated  as  perturbations  on  the  elements  of  the  lunar  orbit 
and  rotation.      Thus,   the  three  components  are  called  libration  in 
longitude    T  ,    in  node    a ,   and  in  inclination    p .      The  libration  in 
longitude  is  a  perturbation  on  the  celestial  longitude  of  the  moon  (can 
be  considered  as  a  correction  to  lunar  time),   libration  in  node  on  the 
celestial  longitude  of  the  node,   and  libration  in  inclination  on  the 
inclination    I    of  the  moon's  equator.      These  perturbations  are  small, 
never  exceeding  2'  of  arc.      They  can  be  expressed  as  follows: 

T       =       S  at  sin  (ng  +  mg'  +  pF  +  qD) 

Icr      =       It  +Sbi  sin(ng  +mg'  +  pF  +qD) 

p       =       Scj  cos  (ng  +  mg'  +  pF  +  qD) 

Where  g  and  g'  are  the  mean  anomaly  of  the  moon  and  sun  respec- 
tively  (dependent  on  time). 

F      =       g  +  w 

D      =       g-g'+w-w1 

Where  w  and  w'  are  the  arguments  of  the  perigrees  of  the  moon  and 
sun  (dependent  on  time),     n,   m,   p,   and  q   and  the  coefficients    a,  b,   and 
c    are  given  in  Table  11. 
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Symbols 

Series 

Argument 
n            m 

Multipliers 

p       q 

Coefficients  in 
Arc  Seconds 

T 

sin 

0 

0 

2 

-2 

ai  = 

+    1.7 

0 

1 

0 

0 

+  91.7 

1 

-1 

0 

-1 

-     1.2 

1 

0 

0 

-2 

+    4.2 

1 

0 

0 

-1 

-     3.5 

1 

0 

0 

0 

-  16.9 

2 

-1 

0 

-2 

+    1.0 

2 

0 

-2 

0 

+  15.3 

2 

0 

0 

-2 

+  10.0 

la 

sin 

0 

0 

2 

-2 

bi  = 

-    3.2 

0 

0 

2 

0 

-  10.6 

1 

0 

-2 

0 

-  23.8 

1 

0 

0 

-2 

+    2.5 

1 

0 

0 

0 

-100.7 

P 

cos 

0 

0 

2 

-2 

Ci   = 

-     3.2 

0 

0 

2 

0 

-  11.0 

1 

0 

-2 

0 

+  23.9 

1 

0 

0 

-2 

-     1.9 

1 

0 

0 

0 

-  98.5 

Table  11.    -   Physical  Libration  Argument 
Multipliers  and  Coefficients 
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APPENDIX  D 
BLACKSHEAR  SPHERICAL  HARMONIC   COEFFICIENTS 

Table  12  below  gives  the    C    and    S    coefficients  as  determined  by 
Blackshear  in  1969.      The  coefficients  are  to  the  13th  order  and  degree. 
These  coefficients  are  connected  with  satellite  dynamics,   with  the 
potential    V    written  in  the  form  [9]  : 

00 

Km  N        a 

V  =   ~  £l  "  S  S     (")"  [JnMRnm  (6,  X)  +KNMSNM  (0,X)]  } 

1  N=l  M=0       L 


N 

M                                 C 

S 

0 

0            \  .000057 35c\\^  +  00 

0. 
0. 

j 

: 

2 

0         -2.0707075166E-04 

3 

0         -6.3026014213E-06 

0. 
0. 
0. 

4 

0             1.9378612685E-05 

5 
6 

0         -7.4593773641E-06 
0         -1.0784697393E-06 

0. 
0. 

6. 

7 

0            2.4033309155E-05 

8 

0         -2.6548898776E-05 

9 

0          -1.5429544143E-06 

0. 

10 

0             5.6344945025E-05 

0. 

11 

0         -2.4602196863E-05 

0. 

12 

0             3.2989171030E-05 

0. 

13 

0         -5.7723614632E-05 

0. 

:06~~ 
-06 

-06 

-06 
-06 
-06 
-06 
-06 
-06 

2 

1         -4.4251194134E-07 

-4.5725447265E- 

3 

1             2.4370024362E-05 

__  2.  3006376707E- 

6.8017  76  5  237E- 

-1.2107811323E- 

2.4894914102E- 

8.0724555222E- 

'  -1.9291~37C2~52E- 

4 

1          -5.7323804405E-06 

5 

1             1.3607744245E-06 

6 

1          -3.6076874485E-07 

7 

1             7.8250002876E-06 

8 

1             3.9152106632E-06 

9 

1             3.1134712059E-06 

-5.33S0336547E- 
-4,  3369907744E- 

10 

1             9.8656663114E-07 

11 

1         -3.01 81405386E-06 

zXa.  1 6  2  2.9 1  39  6  9  E  - 
1.7411943335 E- 

~  -5.75S4594581E- 

-Q5L  . 

-06 

12 

1          -7.4916591513E-06 

13 

1         -1.2121962226E-06 

-06     | 

Table  12.    -  Blackshear  Coefficients 
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6 

6 

-2.0267463326E- 

•09 

3.8994163212E- 

-09 

7 

6 

-1.0459307494E- 

■09 

-3.0269761767E- 

10 

8 

6 

-3.2667897321E- 

10 

-8.8897882759E- 

-10 

9 

6 

1.7328428254E- 

10 

3.2542881147E- 

10 

10 

6 

9.5432215370E- 

11 

6.065CC9C402E- 

-11 

11 

6 

4.3477571496E- 

11 

-1.7188426198E- 
._  ._3 .  7 2.7 Q J  3.4.3  98 E^ 

-11 

.12- 

_6__ 

.  jt*  2  4  828  8607  8_E  - 

ill- 

13 

6 

4.3675063968E- 

■12 

-1.  5336991203E- 

-11 

7 

7 

1.2026816779E- 

-10 

-3.2233006465E- 

-10 

8 

7 

6. .01  88.7  86.28  4  E: 

11 

3.5932703797E- 

-11 

9 

7 

-1.3023550339E- 

-11 

1.0561081783E- 

-10 

10 

7 

-1.6122513843E- 

11 

-2.3325805068E- 

-11 

11 

7 

-3.0064578083E- 

-12 

-1.0710145347E- 

-11 

12 

7 

-1.5484033273E- 

12 

-2.3107720838E- 

-12 

13 

7 

-1.4899022434E- 

-12 

-2.4044711145E- 

-12 

8 

3 

-1.7001116994E- 

11 

.     3.0033533379E- 

11 

9 

8 

-2.8584350275E- 

■12 

7.8834743207E- 

-14 

10 

8 

1.5891710805E- 

■12 

-6.5597619378E- 

-12 

11 

8. 

1.0177719212E- 

-12 

-JU&Lfil  698 197E: 

-.13 

12 

3 

1.0400719186E- 

13 

4.4891103458E- 

-13 

13 

8 

-2.  1168402744E- 

14 

5.5719251933E- 

-14 

9 

9 

9.6176983415E- 

13 

-  1.3039134543E- 

-12 

10 

9 

-1.9696B87071E- 

14 

-2ol938879453E- 

13 

11 

9 

..-  JU.  56  8  2  88  2  50  4  E- 

■13 

3.5220136268E- 

-A3. 

12 

9 

-7.0712057255E- 

14 

-3.  553489C024E- 

.14 

13 

9 

4.6826469874E- 

-16 

-1.9224377776E- 
4.6557201126E- 

14 

10 

10 

-3.6774006812E- 

-14 

-14 

11 

10 

-8.1140071909E- 

-15 
-15 

1.0112926706E- 
"-  1  ."13620893  8  1E: 

14 

12 

10 

4.9625158335E- 

14 

13 

10 

2.  13S0269179E- 

-15 

U4722610146E- 

-15 

11 

11 

3.3647966771 E- 

-16 

-1.1213235338E- 

15 

12 

11 

7.4957632120E- 

-17 
:16 

-2.4793959309E- 

16 

13 

11 

-2.0181618293E- 
"-3~.62746~9163  3E: 

1.7194186339E- 

16 

12 

12 

-18 

•     -1.5526294510E- 

18 

13 

12 

-1.8762164620E- 

-18 

-8. 1189021267E- 

-18 

13  . 

13 

9.5540605576E- 

-19 

-1.7059499010E- 

-19  . 
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APPENDIX  E 

OPTIMUM   POSITIONS  FOR  MINIMUM 
STAR  COORDINATE  ERRORS 

Stars  which  yield  the  minimum  error  in  transformation  were  found 
to  possess  the  coordinates  of: 

Lunar  Declination  =  0 

o  ° 

Lunar  Right  Ascension     =  0    or  180 

Thus  this  point  will  fall  on  the  ecliptic  at  the  position  of  the  nodes. 
Figures  10  and  11    (from  [16] )  give  the  positions  of  these  stars  in  earth 
coordinates  of  right  ascension  and  declination  for  the  time  January  0.5 
1971,    January  0.5  1972,    January  0.5  1973  and  January  0.5   1974.      The 
motion  of  these  points  is  in  direct  relationship  with  the  motion  of  the 
lunar  nodes. 
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